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Grain  size  is  one  of  the  factors  which  influence  mechanical  properties  of  metals  like  strength  and 
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metal  samples  with  identical  grain  size  distribution  exponents,  where  the  shifts  of  the  log-log 
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I.  Introduction 


Grain  size  is  one  of  the  factors  influencing  mechanical  properties  of  metals  like  strength  and 
fracture  toughness,  and  has  been  widely  used  for  monitoring  product  quality  during  manufacturing 
processes.  Ultrasonic  waves  propagating  in  polycrystalline  materials,  such  as  metals,  are  subject 
to  attenuation  often  dominated  by  grain  boundary  scattering.  The  strong  correlation  between 
ultrasonic  attenuation  and  grain  size  has  been  well  established  and  extensively  utilized  in  recent 
years  for  grain  size  estimation  because  it  is  practical,  relatively  inexpensive,  and  most  significantly 
nondestructive.  However,  the  fundamental  mechanisms  underlying  the  empirical  correlation 
between  ultrasonic  scattering  measurements  and  microstructure  are  highly  complex  and  lack  a 
satisfactory  theoretical  model.  Current  practical  techniques  for  material  characterization  and  grain 
size  estimation  are  generally  based  on  simplified  scattering  models  and  often  can  serve  only  limited 
purposes. 

The  importance  of  grain  size  estimation  for  industrial  applications  warrants  the  investigation 
of  alternative  methods  of  relating  ultrasonic  attenuation  to  grain  size.  The  property  known  as 
scaling,  which  can  be  broadly  defined  as  adjustments  according  to  a  proportion,  will  be  defined 
both  mathematically  and  through  examples.  The  application  of  scaling  behavior  for  the 
interpretation  of  ultrasonic  attenuation  experimental  data  will  be  investigated  and  related  to  the 
microstructural  properties  of  the  material  samples.  Although  the  average  grain  size  contains 
valuable  information  related  to  the  mechanical  properties  of  the  material,  the  distribution  of  grain 
size  can  significantly  affect  the  interpretation  of  experimental  results  [Smith  et.  al.  1981].  It 
follows  that  the  knowledge  of  the  grain-size  distribution  can  be  valuable  information  for 
microstructure  characterization.  Relevant  scaling  concepts  will  be  reviewed  and  then  used  to  link 
the  wavelength  dependence  of  the  attenuation  coefficient  extracted  from  experimental  data  directly 
to  the  grain-size  distribution. 

The  contributions  of  this  work  begin  with  the  hypothesis  of  a  new  scaling  relationship 
between  the  grain-size  distribution  and  attenuation.  The  outcome  of  the  theoretical  development  is 
a  straightforward  relationship  between  an  exponent  to  describe  the  grain-size  distribution  and  the 
power-law  dependence  of  attenuation  on  wavelength.  This  establishes  the  first  straightforward 
relationship  between  the  grain-size  distribution  function  and  the  exponent  that  characterizes  the 
dependence  of  attenuation  on  wavelength,  linking  the  microstructural  properties  to  the  dynamical 
response. 

An  experiment  to  verify  this  relationship  for  nickel  samples  is  outlined  in  section  V,  and  the 
results  are  presented  in  section  VI.  A  second  scaling  relationship  is  examined  for  the  annealing 
durations  in  section  VI. D.  The  conclusions  are  discussed  along  with  suggestions  for  future 
directions  in  section  VII. 

The  work  reported  here  is  a  result  of  interdisciplinary  cooperation  of  laboratories  from 
Drexel  University:  the  Signal  Processing  Laboratory,  supervised  by  Drs.  Nihat  Bilgutay  and  Banu 
Onaral,  and  the  Metallographic  Laboratory,  supervised  by  Dr.  Samuel  Nash,  as  well  as 
laboratories  supervised  by  Dr.  E.  R.  Generazio  from  the  NASA  Lewis  Research  Center. 

II.  Review  of  Grain  Size  Considerations 

A.  Measurement  of  Average  Grain  Size  and  Grain  Size  Distributions 

The  accurate  measurement  of  grain  size  is  a  continuing  topic  of  research,  even  when  it  is 
accomplished  destructively  using  polishing,  etching  and  photomicrographs  of  the  grain  boundary 
structure.  The  most  commonly  used  techniques  for  average  grain  size  measurement  are  the  line 
intercept  method  and  the  American  Society  for  Testing  Materials  (ASTM)  comparison  method 
[DeHoff  19®  p.  233],  The  line  intercept  approach  uses  a  grid  placed  over  the  photomicrograph 
and  determines  the  average  number  of  times  test  lines  intersect  the  boundaries.  The  alternative 
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method  compares  the  photomicrograph  to  ASTM  grain  size  charts.  Both  of  these  techniques  are 
highly  subjective  measurements;  line  intercept  results  have  been  shown  to  vary  as  much  as  50  per 
cent  [DeHoff  1968  p.  250]. 

This  work  requires  the  measurement  of  the  grain-size  distribution,  not  just  the  average 
grain  size.  The  distinction  between  average  grain  size  and  grain-size  distribution  measurement 
broaches  an  important  criticism  of  past  work.  Ultrasonic  attenuation  is  proportional  to  the  volume 
of  the  scatterer  for  one  scattering  region,  the  Rayleigh  region.  In  cases  where  there  is  a  distribution 
of  grain  sizes,  the  third  statistical  moment  of  the  grain-size  distribution  should  be  used  in  the 
analysis  of  this  type  of  scattering.  Researchers  have  used  the  third  power  of  the  average  for  the 
third  moment  [Klinman  &  Stephenson  1981],  an  equality  that  may  not  be  valid  when  a  grain-size 
distribution  exists.  In  some  cases,  this  use  of  the  average  may  be  due  to  the  investigators  using  the 
intercept  method  for  grain  size  determination,  which  provides  only  the  average  size  [Mercier  1975, 
Klinman  &  Stephenson  1981,  Saniie  &  Bilgutay  1986]. 

The  effect  of  the  distribution  of  grain  sizes  on  mechanical  properties  of  the  metal  has  been 
suggested,  specifically  that  conventional  experimental  interpretation  may  not  be  correctly 
considering  the  “uncontrolled  and  unrecognized  variations  in  the  grain-size  distributions.”  [Rhines 
&  Patterson  1982]  The  authors  also  state  that  ‘The  near  absence  of  recognition  of  these  features 
and  of  their  effects  on  properties  is  understandable  in  view  of  the  great  amount  of  labor  required  in 
the  experimental  measurement  of  grain  volume  distribution,  together  with  the  fact  that  the 
distributions  are  not  particularly  noticeable  in  microstructures  seen  two-dimensionally.”  [Rhines  & 
Patterson  1982] 

One  recent  innovation  in  the  determination  of  the  grain-size  distribution  is  the  use  of  two- 
dimensional  Fourier  transforms  of  an  image  of  the  grain  boundaries  [Generazio  1988b].  This 
method  is  briefly  reviewed  below  since  it  is  used  for  the  experimental  work  reported  here. 

A  photomicrograph  of  the  grain  boundaries  is  first  digitized  using  a  512  by  512  pixel  array, 
circularly  masked,  and  then  the  grain  boundaries  are  enhanced  digitally.  The  grains  are  tone-pulse- 
encoded  for  Fourier  domain  measurement  of  the  grain  sizes  as  follows:  beginning  at  the  center  of 
the  image  and  proceeding  outwards  along  a  radial  line,  3  periods  of  a  sine  wave  are  fitted  into  the 
width  of  each  grain  crossed  by  the  radial  line.  This  creates  a  tone-pulse-encoded  set,  whose 
Fourier  transform  reveals  the  density  function  for  the  widths  (or  size)  of  the  grains.  A  one¬ 
dimensional  example  of  tone-pulse  encoding  is  shown  below  in  Figure  1.  This  process  is  carried 
out  for  a  fixed  number  of  radial  lines  circling  the  micrograph.  A  two-dimensional  density  function 
dependent  on  size  and  angle  is  created.  If  the  grain-size  distribution  is  isotropic,  then  the 
histogram  will  be  independent  of  angle. 

A  test  case  for  this  method  [Generazio  1988b]  utilized  an  image  with  a  known  distribution 
of  circles.  The  algorithm  correctly  measured  the  size  distribution  for  the  larger  grains,  but  it 
deviated  for  the  smaller  grains.  This  is  easily  understood  by  analyzing  how  the  circles  are  tone- 
pulse  encoded.  Figure  1  is  clearly  an  ideal  case:  the  circles  are  encoded  along  their  diameter,  the 
longest  chord  of  the  circle.  However,  as  the  size  distribution  measurement  is  taken  along  another 
radial  line,  the  chord  will  not  be  the  diameter,  and  the  measurement  will  give  additional,  misleading 
weighting  to  smaller  sized  grains.  Consequently,  the  size  measurement  for  the  smaller  grains  is 
expected  to  have  a  higher  degree  of  uncertainty  than  the  measurement  for  the  larger  grains. 

The  percent  uncertainty  in  the  grain  size  measurement  [Generazio  1988b]  can  be  calculated 
as  follows: 


percent  uncertainty  = 


3  pixels 


*  100  per  cent 


(1) 


where  d  is  the  grain  size  and  3  pixels  is  the  experimental  uncertainty  based  on  the  size  erf'  the  grain 
boundary  on  the  micrograph.  This  expression  provides  another  reason  why  the  measurement  is 
more  accurate  for  larger  grain  sizes;  the  smaller  the  grain,  the  larger  the  percent  uncertainty. 
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Additionally,  since  the  Nyquist  rate  [Oppenheim  &  Willsky  1983,  p.  519]  must  be  satisfied,  the 
smallest  number  of  pixels  that  can  be  used  to  encode  a  three-cycle  sine  wave  is  ten.  This 
observation  sets  the  lower  limit  for  the  grain  size  that  can  be  resolved  using  this  technique. 

II. B.  Distribution  Functions  Used  For  Grain  Size 

It  will  be  shown  that  the  attenuation  coefficient  can  be  related  to  the  grain-size  distribution 
function,  N(D)  defined  as 

N(D)  dD  =  the  number  of  grains  of  size  between  D  and  D  +  dD 

It  is  often  difficult  to  assign  a  mathematical  formula  for  the  distribution  function,  especially  when 
only  a  few  data  points  are  known,  as  is  the  case  for  many  grain  size  measurements. 

The  log-normal  distribution  has  been  used  for  grain  size  in  metals,  based  on  measurements 
with  twelve  values  [Feltham  1957],  and  [Seemann  &  Bentz  1954]  (German  text)  as  reported  in 
[Papadakis  1961].  It  is  interesting  to  note  that  one  author  [DeHoff  1965]  indicates  that  based  on 
experimental  results  the  distribution  must  be  skewed  and  that  the  log-normal  is  not  the  only 

possible  distribution  function;  it  simply  “provides  a  versatile  functional  form _ In  the  absence 

of  additional  information,  it  is  recommended  that  this  form  be  used  . . . .”  [DeHoff  1965]. 

The  log-normal  distribution  has  been  used  for  many  other  applications,  including  body 
weights,  rainfall,  and  electrical  resistance  in  materials  with  random  scatterers  [West  &  Shlesinger 
1989].  It  is  one  of  a  class  of  distributions  used  for  systems  comprised  of  many  components.  It 
has  been  shown  that  there  is  a  connection  between  the  log-normal  and  another  distribution 
applicable  to  many-component  systems,  the  power-law  distribution.  The  upper  tail  end  of  a  log¬ 
normal  distribution  (the  last  two  to  three  percentile  region)  can  be  approximated  by  an  inverse 
power-law  [West  &  Shlesinger  1989],  as  can  a  distribution  with  a  large  variance  [Montroll  and 
Shlesinger  1983].  While  an  inverse  power-law  expression  was  used  for  the  grain-size  distribution 
function  in  this  work,  it  is  interesting  to  keep  in  mind  these  relationships  between  the  log-normal 
and  power-law  functions. 

There  is  yet  another  connection  between  the  log-normal  and  power-law  distributions.  The 
cumulative  frequency  of  the  grain  weight,  i.e.,  the  percentage  of  grains  which  are  less  than  or 
equal  to  a  given  weight,  is  often  represented  on  the  log-log  scale  [Rhines  &  Patterson  1982, 
Doherty  1984].  The  linearity  of  this  figure  is  designed  to  show  that  the  volume  distribution  is  log¬ 
normal;  perhaps  a  simpler  way  of  interpretation  is  to  represent  the  volume  distribution  with  a 
power  law. 

III.  Ultrasonic  Attenuation 
A.  Introduction 

The  ultrasonic  wave  traveling  through  a  solid  is  subject  to  scattering  and  absorption  which 
results  in  energy  losses.  As  the  ultrasonic  burst  advances  through  homogeneous  material,  its 
amplitude  is  attenuated  as 

P  =  Poe"ad  (2) 

where  p0  is  the  initial  pressure,  p  is  the  pressure  at  d,  the  distance  traveled,  and  a  is  the  pressure 
attenuation  coefficient  in  Nepers  per  centimeter  (Np/cm)  [Krautkramer  &  Krautkramer  1983  p.  . 

106].  In  general,  the  attenuation  coefficient  can  be  a  function  of  wavelength  (X)  and  size  of  the 
scatterer  (D)  and  can  be  expressed  as 

a  =  aa+aK  (3) 
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where  is  the  absorption  term  and  is  the  scattering  loss  associated  with  grain  boundaries. 
Since  wavelength  and  frequency  (0  are  related  through  the  velocity  of  propagation  c  (cm/sec): 

c  =  f  X  (4) 

the  attenuation  can  also  be  expressed  in  terms  of  its  frequency  dependence. 

Absorption  in  metals  is  caused  primarily  by  the  inelastic  behavior  of  the  materials, 
dislocation  damping,  and  magnetic  domain  damping  [Papadakis  &  Reed  1961,  Smith  1987],  The 
dependence  of  the  absorption  loss  on  frequency  is  usually  presumed  to  be  linear  [Roth  1948]. 
Many  researchers  have  assumed  that  the  grain  scattering  loss  at  ultrasonic  frequencies  is  so  large 
relative  to  the  absorption  loss  that  the  latter  is  negligible  [Papadakis  1970,  Serabian  &  Williams 
1978,  Bozorg-Grayeli  1981  p.  16,  Smith  et.  al.  1981].  One  attempt  to  justify  this  assumption 
experimentally  [Wegel  &  Walther  1935]  used  different  materials,  most  notably  copper.  Ultrasonic 
attenuation  measurements  were  made  before  and  after  annealing.  Since  the  copper  before 
annealing  has  little  or  no  grain  boundaries,  the  dissipation  of  energy  was  due  to  absorption.  After 
annealing,  the  dissipation  immediately  increases  dramatically  (by  a  factor  of  more  than  16)  due  to 
recrystallization  and  the  formation  of  grain  boundaries.  Unfortunately,  there  is  a  drawback  to  the 
use  of  this  experiment  as  a  means  of  justifying  the  assumption  that  grain  boundary  scattering 
dominates  ultrasonic  attenuation.  The  dislocations  in  the  crystalline  structure  also  change  during 
annealing,  and  they  contribute  to  the  absorption  term  of  the  attenuation.  Thus,  it  is  impossible  to 
separate  the  increase  of  the  attenuation  after  annealing  due  to  either  formation  of  grain  boundaries 
or  change  of  the  dislocations.  This  work  is  therefore  left  with  the  assumption  that  the  attenuation  in 
the  metal  samples  is  due  to  grain  boundary  scattering  alone;  thus  the  notation  for  the  attenuation 
will  be  a(X,D)  [i.e.,  a(X,D)  »  a^X.D)]. 

A  term  linearly  dependent  on  frequency  has  been  previously  added  to  the  attenuation 
expression  for  some  applications  [Mason  &  McSkimin  1948].  This  term  assumably  accounts  for 
the  power  dependence  deviating  from  four,  as  expected  for  Rayleigh  scattering.  In  other 
developments,  an  f2  term  was  used,  again  based  on  the  frequency  dependence  differing  from  f4 
[Kumar  &  Govil  1980], 

III.B.  Attenuation  Theory 

1.  Single  Grain  Size 

Ultrasonic  scattering  formulas  have  been  studied  and  classified  into  distinct  scattering 
regions  according  to  the  ratio  of  the  sound  wavelength  to  the  size  of  the  grains.  This  section 
reviews  some  important  results  in  detail  for  material  with  a  single  grain  size  or  an  average  grain 
size.  The  review  of  the  extension  of  the  theory  to  a  grain-size  distribution  [Smith  1982]  follows  in 
Section  III.B.2.  The  assumptions  used  in  the  derivations  can  include: 

•  grains  are  weakly  anisotropic  (dependence  on  direction) 

•  there  is  no  preferred  orientation  for  grains 

•  grains  are  equiaxed  or  spherical 

•  material  contains  single-crystal  grains  of  a  single  phase  with  no  voids  or 
inclusions 

The  theoretical  development  based  on  these  assumptions  relates  the  attenuation  due  to  grain 
boundary  scattering  to  wavelength,  and  in  many  cases  to  the  dependence  on  elastic  constants  and 
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sound  velocities  of  the  material.  It  is  reviewed  in  the  section  below  that  the  attenuation  is  expected 
to  vary  as  a  function  of  wavelength  (X)  and  scatterer  size  (D)  as  follows: 


Rayleigh  Region: 

a(X,D)  =  Cr  D3  X4, 

X»D 

Stochastic  Region: 

a(X,D)  =  Q  DX'2, 

X  ~  D 

Diffusive  Region: 

a(X,D)  =  Cd/D, 

X«D 

(5) 

where  Cp  Cs,  and  Cj  are  constants  of  the  material.  For  spherical  grains,  D  is  the  diameter;  in  the 
case  of  nonspherically  shaped  grains,  D  is  interpreted  to  be  the  largest  measurement  of  the  grain. 

There  have  been  numerous  studies  of  the  theoretical  developments  resulting  in  the 
expressions  for  the  attenuation  in  Equation  5.  Selected  references  are  cited  below,  providing 
examples  of  approaches  and  assumptions  used  to  derive  the  attenuation  for  the  different  scattering 
regions: 


Table  1:  References  for  Theoretical  Approaches  to  Attenuation 


a)  X  »  D,  single  scatterer,  neglecting  mode 
conversion 

Rayleigh  1894 

b)  X  »  D,  multiple  scatterers,  neglecting 
mode  conversion 

Mason  &  McSkimin  1947 

c)  X  »  B,  multiple  scatterers,  using  mode 
conversion 

Lifshitz  &  Parkhomovskii  1950, 
Merkulov  1956,  Bhatia  1959a, 
1967,  Bhatia  ft  Moore  1959b 

d)X«  B,  multiple  scatterers 

Mason  &  McSkimin  1948 

e)  all  relative  values  of  X  and  D 

Roney  1950 

0  statistical  signal  processing 

Saniie  &  Bilgutay  1986,  Saniie 
et  al.,  1988 

III.B.2.  Grain-Size  Distributions 

The  effect  of  neglecting  the  grain-size  distribution  can  be  seen  in  the  interpretation  of 
experimental  results.  Some  investigators  have  postulated  that  the  ultrasonic  attenuation  is 
dominated  by  Rayleigh  scattering  and  have  used  a  fourth  power  dependence  of  attenuation  on 
frequency  in  order  to  find  the  average  grain  size  [Kumar  ft  Govil  1980,  Mercier  1975,  Papadakis 
1964,  1965a,  1965b,  Roderick  ft  Truell  1952,  Smith  et.  al.  1981],  despite  measurements 
demonstrating  that  the  power  dependence  has  a  noninteger  exponent  between  2  and  4  [Mason  ft 
McSkimin  1947,  Papadakis  1970].  In  other  words,  researchers  have  used  the  following 

relationship  for  the  average  grain  size  (B)  and  frequency  (0 

a(f)~(B)3f4  (6) 

with  a  certain  degree  of  success,  but  this  approach  has  proved  unsatisfactory  for  some  samples  for 
two  possible  reasons.  First,  using  the  third  power  of  the  average  grain  size  may  not  be  valid  when 
there  is  a  distribution  expression  for  the  grain  sizes.  When  a  grain-size  distribution  exists,  the 
attenuation  is  a  function  of  the  third  statistical  moment  Second,  Rayleigh  scattering  requires  that 
the  wavelength  is  much  greater  than  the  size  of  the  scatterers.  Although  the  wavelengths  chosen 
may  be  greater  than  the  average  grain  size,  they  may  in  fact  be  comparable  in  size  to  the  largest 
grain  size,  causing  other  scattering  processes  to  occur.  Thus  one  or  both  of  the  two  terms  of 
Equation  5  can  be  inappropriate  for  the  dependence  of  attenuation  on  grain  size  and  frequency. 


6 


The  most  developed  application  of  grain-size  distributions  to  attenuation  was  given  in 
[Smith  1982],  for  which  both  Rayleigh  and  stochastic  scattering  was  assumed  to  occur.  The 
author  began  with  Roney's  expression  for  scattering  for  all  wavelengths  [Roney  1950  p.  14] 


a  = 


Dpn  K 


(f)  l(2m + l)s,nl6" 


(7) 


and  used  a  probability  density  function  for  the  grain  diameters.  The  total  attenuation  per  unit  length 
due  to  all  scatterers  was  found  neglecting  multiple  scattering,  and  then  the  author  then  used  a  log¬ 
normal  grain-size  distribution. 


IV.  Scaling  Approach  to  Wavelength  and  Grain  Size  Dependence  of  Ultrasonic 
Attenuation 


A.  Introduction  to  Scaling  Approach 

The  objective  of  this  work  is  to  introduce  and  apply  a  novel  approach  to  ultrasonic 
attenuation  based  on  scaling  analysis.  Scaling  can  be  defined  as  a  property  describing  objects  or 
systems  that  are  invariant  under  certain  transformations  of  scale  [Mandelbrot  1983  p.  18].  The 
discussions  in  this  section  demonstrate  that  power-laws  are  one  appropriate  expression  for  scaling 
properties.  Many  research  topics  involve  the  study  of  both  the  structural  and  dynamical  scaling 
properties  of  a  system;  the  research  presented  here  falls  into  this  class  of  investigation.  Fractal 
concepts  are  closely  associated  with  scaling  because  the  irregularity  or  fragmentation  for  many 
fractal  objects  are  identical  at  all  scales. 

The  idea  for  the  scaling  work  originated  with  the  experimental  observation  that  the 
wavelength  dependence  of  ultrasound  attenuation  is  an  inverse  power-law,  with  a  noninteger 
exponent  [Mason  &  McSkimin  1947,  Papadakis  1970,  Kumar  &  Govil  1980].  This  is  in  contrast 
to  the  classical  derivations  in  Equation  5,  which  demand  a  dependence  on  an  integer  (0,  2  or  4) 
power  of  the  wavelength.  Scaling  concepts  provide  an  approach  to  the  interpretation  of  a 
measurement  that  exhibits  this  noninteger  power  dependence,  and  can  result  in  the  ability  to 
characterize  structures  or  dynamics  in  terms  of  straightforward  scaling  measures  that  can  yield 
models  with  a  small  number  of  parameters. 

Self -similarity  is  a  special  kind  of  scaling,  when  a  piece  of  an  object  is  geometrically  similar 
to  the  whole  [Mandelbrot  1983  p.  18].  One  consequence  of  self-similarity  is  that  the  properties  of 
a  self-similar  geometric  object  are  indistinguishable  as  a  function  of  resolution,  also  known  as 
length  scale  [Or bach  1986].  No  matter  how  much  an  image  of  the  object  is  magnified  or  reduced, 
the  geometric  properties  remain  the  same.  This  condition  is  termed  ‘no  characteristic  scale.’ 

Self-similarity  may  be  defined  in  an  exact  or  statistical  sense.  The  definition  for  self¬ 
similarity  above,  namely  the  property  that  a  piece  erf  an  object  is  geometrically  similar  to  the  whole, 
is  exact  self-similarity.  Some  applications  used  the  less  restrictive  definition  of  statistically  self¬ 
similarity:  the  statistical  aspects  are  scale  invariant  [Sander  1987,  Feder  1988  p.  185,  Kaye  1989  p. 
35].  This  work  concentrates  on  statistical  self-similarity.  It  has  also  been  found  that  some 
structures  in  nature  may  be  self-similar  over  a  finite  range  of  length  scales  [Orbach  1986]. 

Concepts  from  scaling  theory  can  be  used  to  justify  the  consideration  of  an  inverse  power- 
law  for  the  distribution  of  grain  sizes.  It  is  reviewed  below  that  if  the  size  distribution  is 
statistically  self-similar,  it  can  be  described  by  an  inverse  power-law,  applicable  when  there  is  no 
characteristic  scale  for  the  measurement  of  the  size  distribution.  To  analyze  the  self-similarity,  the 
scaling  argument  is  applied  to  two  images  of  grains  at  different  magnifications  with  finite 
resolutions:  magnification  1,  and  a  higher  magnification  2.  Figure  2  is  provided  for  illustration 
purposes: 
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Histograms  can  be  made  at  each  magnification  using  Ni(d)  and  N2(d),  defined  as  follows: 

N i(d)  dd  =  number  of  circles  of  diameter  between  d  and  dd  at  magnification  1 

N2(d)  dd  =  number  of  circles  of  diameter  between  d  and  dd  at  magnification  2 

The  question  is  this:  how  does  the  histogram  vary  as  the  magnification  changes?  Unlike  the 
coastline  example,  a  function  rather  than  a  single  value  is  generated  at  each  magnification,  and  an 
appropriate  test  for  possible  scaling  relationships  must  be  chosen.  The  distribution  functions  are 
compared  at  relative  positions  using  the  following  equation: 


N  i(d)  ,  N2(d’) 

RI(37E)  =?W7E) 


(8) 


One  can  think  of  d  as  the  yardstick,  and  ‘b’  as  a  parameter  to  make  the  yardstick  smaller  (if  b  >  1) 
or  larger  (if  b  <  1).  If  the  equality  above  holds,  then  is  dependent  on  the  parameter  ‘b’ 


alone,  not  upon  the  magnification  used.  In  this  example,  the  circles  are  termed  statistically  self- 
similarly  distributed,  and  there  is  no  characteristic  sense  of  scale  for  the  grain-size  distribution. 
The  inverse  power-law  is  the  expression  of  choice  to  represent  such  self-similar  distributions: 
N(d)  ~  d~P.  This  type  of  distribution  has  been  observed  and  discussed  in  light  and  x-ray  scattering 
experiments  for  silica  gels  [Martin  1986,  Martin  &  Wilcoxon  1988]  and  for  the  pore  size 
distribution  in  lignite  coal  [Schmidt  1982]. 


IV.B.  Attenuation  Coefficient 

This  section  presents  a  scaling  relationship  between  the  ultrasonic  attenuation  coefficient 
and  the  grain-size  distribution.  As  reviewed  in  Section  III,  attenuation  depends  on  the  wavelength 
(X)  and  the  size  of  the  scatterer;  the  three  scattering  regions  are: 


Rayleigh: 

X»D 

_  Cr  D3 

Stochastic: 

X  =  D 

C„  D 

Diffusive: 

X«D 

«XD<X’D)  =  ^ 

(9) 

where  D  is  the  diameter  or  largest  length  measurement  of  the  grain,  and  Cr,  Cs,  and  Cd  are 
constants  of  the  material.  The  notation  a^D  is  used  in  this  section  to  underscore  the  dependence 
of  attenuation  on  grain  size  and  wavelength.  In  Equation  11,  the  attenuation  is  a  function  of 
wavelength  alone,  and  the  symbol  is  used  to  reflect  this. 

The  theoretical  development  begins  by  defining  new  functions,  called  ‘scaled  attenuations’ 
by  multiplying  the  attenuation  by  the  wavelength  [Generazio  1986]: 


?*' 

^Xd^’^  =  =  Cj.  ( 

=  ®XD^’^  =  f-g 

D  , 
X*1 

X°xD(k,D)  =  a^pfX.D)  =  Cd  (  j-) 

(10) 

The  scaled  attenuations  for  these  three  scattering  regions  are  dependent  on  the  ratio  of  grain 
size  to  wavelength.  The  theoretical  development  has  been  given  which  results  in  a  single 
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expression  for  the  attenuation  for  all  wavelengths  using  the  assumption  of  spherical  grains  [Roney 
1950  p.  14].  Also,  the  expressions  for  attenuation  in  Equation  5  also  depend  on  the  ratio  of  grain 
size  to  wavelength  for  the  three  scattering  regions.  It  is  therefore  assumed  that  the  scaled 
attenuation  is  always  dependent  on  the  ratio  of  grain  size  to  wavelength,  for  all  ratios  of  grain  size 
to  wavelength. 

To  account  for  the  distribution  of  grain  sizes,  the  function  N(D)  is  used,  defined  as 

N(D)  dD  =  the  number  of  grains  of  size  between  D  and  D  +  dD 

Assuming  that  multiple  scattering  can  be  ignored,  the  scaled  attenuation  resulting  from  the 
distribution  of  grain  sizes  is  expressed  as  follows: 


ax(X)  = 


D)  dD 


(11) 


The  sizes  are  assumed  to  be  distributed  with  an  inverse  power-law  with  scaling  exponent  y: 

N(D)  =  KD“y,  0  <  D  <  oo 
00 

a x(k)  =  JKD"Y  a(X,D)  dD  (12) 

One  reason  for  this  assumption  can  be  based  on  statistical  self-similarity  of  the  size 
distribution.  Also,  as  discussed  in  section  II,  the  cumulative  frequency  erf  the  grain  weight  is  often 
represented  linearly  on  the  log-log  scale  [Rhines  &  Patterson  1982,  Doherty  1984],  suggesting 
power-laws  may  be  appropriate  models  for  grain  sizes.  Lastly,  it  was  pointed  out  in  the  same 
section  that  the  log-normal  distribution  can  be  approximated  by  an  inverse  power-law  for 
distributions  with  large  variances. 

The  scaled  attenuation  functions  in  Equation  10  are  dependent  on  the  ratio  This  fact  will 

be  exploited  below,  using  the  notation  a^x)  where  x  =  The  change  of  variables  x  =^-  in 
Equation  12  reveals 

00  00 

Kx"y  K"y  ax(x)  Xdx  =  K 

Since  the  above  inverse  power-law  relationship  between  the  grain-size  distribution  and 
attenuation  has  not  been  previously  examined,  the  convergence  of  the  integral  in  Equation  13  is 

given  below.  This  results  in  bounds  for  the  exponent  y.  The  integral  is  expressed  using  different 
intervals  of  x: 


/' 


x  Y  ax(x)  dx 


(13) 


ax(X)  =  K  k~(r 


a„(x)  dx  + 


*2 


ax(x)  dx  + 


ax(x)  dx 


where  *2  >  xi  *i  *  1  and  *2  *  L 


(14) 
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The  first  integral  in  Equation  14  tests  the  convergence  for  small  x.  Observing  from 
Equation  10  that 

3 

a^x)  ~  x  ,  x  «  1 

the  first  integral  in  Equation  14  is 


xi 

Xj  «  1  J\~ Y  x3  dx  = 

which  converges  if  4  >  y.  Conversely,  for  large  x 
a^x)  ~  x'1,  x  »  1 

and  the  third  integral  in  Equation  14  becomes 


o 


*2  *  1 


(15) 


(16) 


which  converges  when  y  >  0.  Therefore,  the  bounds  for  the  scaling  exponent  are  0  <  y  <  4  in 
order  for  the  integral  in  Equation  13  to  converge.  The  wavelength  dependence  of  a^A)  is 


aJt(X)  ~  A"^Y  0  <  y  <  4 

and  the  expression  for  the  attenuation  coefficient  becomes 


(17) 


=  A,'1  a^A)  =  K  A  Y  Jx  Y  a,(x)  dx  ,  0  <  y  <  4  (18) 

Therefore  the  exponent  for  wavelength  dependence  of  attenuation  is  determined  by  the 
exponent  for  the  grain-size  distribution.  In  other  words, 

the  grain-size  distribution  is  the  dependence  of  attenuation  on 

an  inverse  power  law  with  =*■  wavelength  is  an  inverse  power 

exponent  y  law  with  exponent  y 


It  has  been  proven  that  the  reverse  is  also  true:  if  the  dependence  of  attenuation  on 
wavelength  is  an  inverse  power-law  with  exponent  y  then  the  grain-size  distribution  is  also  an 
inverse  power-law  with  exponent  y  [Nicoletti  1991]. 

In  summary,  the  work  presented  in  this  section  proposes  a  power-law  relationship  between 
the  proposed  grain-size  distribution  function  and  the  resulting  functional  dependence  of  attenuation 
on  wavelength.  Additionally,  this  relationship  is  shown  to  be  valid  if  the  development  starts  with 
tl  ;  assumption  of  an  inverse  power-law  for  the  grain-size  distribution  or  the  attenuation.  This 
establishes  the  first  straightforward  relationship  between  the  grain-size  distribution  function  and  the 
exponent  which  characterizes  the  dependence  of  attenuation  on  wavelength,  linking  the 
microstructural  properties  to  the  dynamical  response. 

IV.C.  Limits  of  Inverse  Power-Law  Distribution 

This  section  analyzes  the  practical  limitations  of  die  results  shown  in  section  V.B, 
demonstrating  the  effect  of  the  choice  of  wavelength  on  the  measured  attenuation  exponent.  By 
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using  the  assumption  that  an  inverse  power-law  distribution  is  retained  within  a  limited  region 
(Di  <  D  <  D2),  the  expression  for  the  scaled  attenuation  becomes: 


Di  D2  00 

a^X)  =  D)  a^D.X)  dD  +jfKD~Y  dD  +  ^(D)  dD 


(19) 

If  the  number  of  grains  with  size  less  than  Di  or  greater  than  D2  is  negligible,  then  the  first  and 
third  integrals  can  be  assumed  to  be  insignificant  in  comparison  to  the  second  integral  and  the 
scaled  attenuation  can  be  written  as 


D2 

a^X)  =  JKD"Y  dD  (2°) 

and  after  rewriting  in  terms  of  x  =  — , 


D2/X  D2/X 

a^X)  =  |kx'Y  X~Y  a^x)  Xdx  =  K  X“*Y  fx-Y  ax(x)  dx  (21) 

Dj/X  Di Tk 

Three  cases  for  the  ratio  of  wavelength  to  Dj  and  D2  are  analyzed  below: 

Case  1:  X  «  Di  (Diffusive  scattering  rally) 

Case  2:  X  »  D2  (Rayleigh  scattering  rally) 

Case  3:  Di «  X  «  D2  (All  three  scattering  modes  simultaneously) 


Case  1  (X  «  Dj) 

If  X  «  Di  then  the  attenuation  should  be  due  to  diffusive  scattering  alone  and  the  exponen  t 
for  the  wavelength  dependence  of  the  resulting  attenuation  should  be  zero  (independent  cr 
wavelength).  From  Equation  10 

x  »1  Vx)  =  ^  (22) 


D 


where  x  =  ~  and  the  expression  in  Equation  21  becomes 


D2/X 


ax(X)  =  K  X"(Y  _1)  fx 

Dil 


-Y  Ch 


‘YTdx 


-Y 


~  X 


Thus  the  attenuation  is  independent  of  wavelength,  as  expected: 
a^X)  ~  X°,  X  «  Di 


(23) 

(24) 
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Case  2  (X  »  D2) 

If  X  »  D2,  then  the  classical  Rayleigh  limit  should  result.  In  this  instance,  the  expression 
for  the  scaled  attenuation  function  reduces  to  the  case  in  Equation  10  for  small  x: 

x «  1  a,(x)  =  Cf  x3  (25) 

resulting  in  the  following  dependence  of  scaled  attenuation  on  wavelength: 

D2/x 

a^X)  =  KX^'0  fx~Y  x3  dx 
D j/x 


V  /D24~y  -  Pi4'Y\  K(D24'y-  D,^) 

4-Y  1  X4'Y  )  (4-y)X3 

Thus  the  exponent  for  the  dependence  erf  attenuation  on  wavelength  is  four,  as  expected: 

ax(X)  ~  X4,  X  »  D2  (27) 


Case  3  (Di «  X  «  P2) 

When  the  range  for  sizes  is  very  wide  (D2  »  Di),  wavelengths  can  be  carefully  selected 
such  that  Di  «  X  «  D2.  For  this  range  of  wavelengths  all  three  types  of  scattering  (Rayleigh, 

stochastic  and  diffusive)  will  occur.  Also,  Equation  21  is  consistent  with  Equation  13  if 

A» 


approaches  0  and  ^-approaches  infinity.  For  this  case,  the  resulting  dependence  of  attenuation 


coefficient  on  wavelength  is  determined  by  the  exponent  of  the  inverse  power-law  expression  for 
the  grain-size  distribution,  as  in  section  IV.B: 


a^fX)  =  X  *  ax(X)  ~  X  Y,  Di  «  X  «  D2 


(28) 


Therefore,  the  choice  of  wavelengths  for  the  ultrasonic  experiment  is  critically  dependent 
on  the  actual  range  of  grain  sizes  in  the  metal  sample.  If  X  «  Di  (Case  1),  the  measurement  of  the 
attenuation  should  be  independent  of  wavelength  (exponent  =  0).  Conversely  if  X  »  D2  (Case  2), 
then  the  measured  attenuation  exponent  is  expected  to  be  four.  When  Di «  X  «  D2  (Case  3),  the 
measured  attenuation  exponent  is  determined  by  the  grain-size  distribution  exponent  These 
transitions  can  be  viewed  as  crossovers  from  integer  to  noninteger  exponents.  The  objective  of 
this  research  effort  is  the  determination  of  the  grain-size  distribution  exponent,  therefore  the 
wavelengths  used  must  be  Di «  X  «  D2. 

While  this  work  concentrates  on  nondestructive^  measuring  the  grain  size  exponent,  the 
average  grain  size  is  desired  for  many  applications.  If  the  limits  of  the  grain-size  distribution  as 
well  as  the  grain  size  exponent  are  known,  then  any  statistical  moment  can  be  calculated,  including 
the  average  grain  size: 


D  =  K-=-2 


D22'y-D,2’y 


(29) 
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IV.D.  Self-similar  Area  from  Self-similar  Volume  Distributions 

The  grain  size  measurement  is  conventionally  made  from  a  micrograph  of  the  grain 
structure,  an  inherently  two-dimensional  measurement.  The  notation  for  the  distribution  of  sizes 
on  the  cutting  plane  is  Na(d),  where 

Na(3)  dd  =  number  of  circles  of  size  between  d  and  d  +  dd 

The  relationship  between  this  two-dimensional  measurement  and  the  size  of  the  grain  (in  three- 
dimensions)  will  be  explored  in  the  following  section.  If  one  considers  spherical  grains*,  the 
relationship  between  Na(a)  and  the  grain-size  distribution  in  three  dimensions,  N(D)  is  [DeHoff 
1968  p.  155] 


NaO)  =  a 


N(D) 


\Id2  -  a2 


dD 


(30) 


The  expression  in  Equation  30  is  known  as  Abel’s  Equation  [DeHoff  1968  p.  155]. 

If  N(D)  is  statistically  self-similar  then  the  scaling  properties  provide  the  following 
equalities: 


N(D)  =  KD"y  N(bD)  =  a  N(D) 

Y  =  -{^  N(bD)  =  b'YN(D)  (31) 

If  it  is  assumed  that  N(D)  is  statistically  self-similar  then  two  questions  of  interest  arise:  is  N^d) 
also  self-similar?  And  if  so,  with  what  exponent?  To  answer  these  questions,  the  following 
renormalization  transformation  is  performed: 

N^kja)  l  l^N^a)  (32) 

for  kj,  k2  >  0.  Equation  30  is  used,  and  then  a  change  of  variables  (D’  =  D/kj)  is  applied: 


00 


00 


Na(k1d)  =  k 


f  N(P’kt) 

1  J  yjg2D’2  -  kj2a2 


kj  dD’  =  k 


N(D’kt)  dD' 
fP'2-a2 


From  Equation  31, 

NfD’kj)  =  kj^D’) 
and  Equation  33  becomes 


(33) 


(34) 


00 


H^kjd)  =  k^  j 


k^NKD’)  dD' 

VD^-d2 


JN(P’)  dD' 

VD^-a2 


=  kf^'^N^a) 


(35) 


*  This  assumption  has  been  used  for  nickel  grains  (Nicholas  1955] 
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Therefore  the  exponent  for  Ng(d)  is  (y  - 1)  and  Na(d)  can  be  written  as 

Na(d)  =  K’a'(Y'1)  (36) 

This  section  thus  proves  that  if  the  grains  are  self-similarly  distributed,  then  the  circles  seen 
on  micrograph  images  of  the  grain  boundaries  must  also  be  self-similarly  distributed  with  an 
exponent  one  less  than  the  exponent  for  the  volume  distribution.  This  is  the  first  development  of 
its  kind,  demonstrating  the  ease  of  transition  of  between  two  and  three  dimensions  for  power-law 
distributions. 

As  noted  in  section  II,  the  log-normal  distribution  has  been  the  only  function  used  by 
previous  researchers  to  describe  the  grain-size  distribution  N(D).  The  distribution  of  sizes  on  the 
cutting  plane  from  a  log-normal  distribution  of  sizes  using  Abel’s  Equation  is 


NaO)  =  d 


CD 

s 


N, 


tot _  tt-(lnD-Dm)2/2o2 


W2™! 


Vd2  -  02 


dD 


(37) 


where  Dm  =  E[lnD]  and  cr2  =  E[lnD2]-(E[lnD])2.  This  expression  for  the  two-dimensional  size 
distribution  is  difficult  to  evaluate;  furthermore,  it  is  clearly  not  log-normal.  Previous  work  on  the 
relationship  between  the  grain-size  distribution  and  the  attenuation  neglected  this  fact  [Smith  1982]; 
the  two-dimensional  size  distribution  obtained  from  the  micrograph  was  fitted  with  a  log-normal 
curve.  Therefore,  the  assumption  erf  an  inverse  power-law  to  describe  the  three-dimensional  grain- 
size  distribution  has  two  important  consequences.  As  developed  in  section  V.B,  this  assumption 
results  in  a  simple  relationship  for  the  wavelength  dependence  of  attenuation.  Furthermore,  this 
section  demonstrates  that  a  second  outcome  erf  this  assumption  is  a  simple  relationship  between  the 
two-  and  three-dimensional  size  distributions. 


V.  Experimental  Design 
A.  Overview 

This  section  provides  a  brief  overview  of  the  experiments  designed  to  verify  the  scaling 
relationship  between  the  grain-size  distribution  and  ultrasonic  attenuation  as  proposed  in  section  V. 
Two  sets  of  measurements  from  the  same  metal  sample  must  be  made:  1)  the  grain-size  distribution 
from  the  photomicrographs  of  the  grain  boundaries,  and  2)  ultrasonic  attenuations.  The  sample 
preparation  and  attenuation  measurement  set-up  take  into  account  the  practical  constraints  erf  actual 
grain-size  distributions  in  metal  and  the  desire  to  obtain  precise  attenuation  measurement 

This  set  of  experiments  constitutes  the  first  attempt  to  establish  the  proposed  scaling 
relationship  between  the  grain-size  distribution  and  attenuation;  therefore,  much  consideration  was 
given  to  the  use  of  the  simplest  material  and  most  precise  attenuation  measurement  possible. 
Preliminary  ultrasonic  experiments  using  stainless  steel  samples  were  completed  in  the  Signal 
Processing  Laboratory  at  Drexel  University.  The  results  proved  to  be  unsatisfactory  due  to 
diffraction  effects  dominating  the  attenuation  measurement  The  decision  was  made  to  use  higher 
frequencies  (>  10  MHz)  for  the  ultrasonic  measurement  and  to  carefully  prepare  metal  samples  to 
diminish  the  diffraction  as  well  as  reflection  coefficient  effects.  The  attenuation  and  grain-size 
distribution  measurements  were  made  at  the  Lewis  Research  Center  of  the  National  Aeronautics 
and  Space  Administration  (NASA),  under  the  supervision  of  Dr.  Edward  Generazio.  The  metal 
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samples  and  photomicrographs  were  sent  to  the  Lewis  Research  Center,  followed  by  a  visit  by  D. 
Nicoletti  to  the  facility. 

A  pure  metal  rather  than  an  alloy  was  chosen  in  order  to  assure  that  only  grain  boundary 
scattering  is  occurring.  In  alloys,  such  as  stainless  steel,  inclusions  can  form  which  will  scatter  the 
ultrasonic  wave  differently  than  the  grains.  Nickel  was  chosen  because  of  its  inexpensiveness  and 
ready  availability  in  commercially  pure  form  (99.6  per  cent  pure). 

The  metal  processing  steps  include  rolling  and  annealing.  To  satisfy  the  arguments  in 
section  V.C,  these  processing  steps  should  result  in  a  grain-size  distribution  that  is  as  wide  as 
possible.  The  size  distribution  should  also  be  spatially  homogeneous,  and  any  anisotropy  in  the 
shape  of  the  grains  should  be  minimized. 

The  grain-size  distribution  was  measured  using  photomicrographs  of  the  grain  boundaries, 
as  discussed  in  section  II.  Therefore  one  step  of  the  sample  preparation  was  the  necessary 
polishing  and  etching  of  the  surface  to  emphasize  the  grain  boundaries,  which  were  recorded  using 
photomicrographs.  A  final  polishing  prepared  the  metal  surface  for  the  precision  attenuation 
measurement. 

The  attenuation  measurement  was  designed  to  balance  the  constraints  of  diffraction, 
significant  for  large  wavelengths,  and  large  attenuation  (>  5  Np/cm),  significant  for  small 
wavelengths.  Large  attenuation  causes  the  energy  of  the  echoes  to  be  scattered  to  the  extent  that 
their  measurement  is  overpowered  by  the  noise  level,  and  is  significant  for  small  wavelengths 
because  attenuation  is  inversely  proportional  to  wavelength.  Fbr  some  wavelength  ranges  and 
material  preparations,  the  reflection  coefficient  is  wavelength  dependent.  Rather  than  assume 
independence  of  wavelength,  the  reflection  coefficient  was  measured  and  corrected  for  to  eliminate 
its  contribution  to  the  attenuation  measurement 

V.B.  Sample  Preparation 

This  section  details  the  material  processing  steps  as  well  as  the  photomicrographic 
techniques.  The  final  result  is  three  nickel  samples  with  different  annealing  durations.  Nine 
photomicrographs  of  the  three  samples  form  the  data  set  used  for  the  grain-size  distribution 
measurement  (twenty-seven  micrographs  in  total). 

Nickel  200  samples  were  obtained  from  Inco  Alloys  International,  Huntington,  West 
Virginia.  Nickel  200  is  commercially  pure  (99.6%)  wrought  nickel,  with  the  following  maximum 
percentage  chemical  composition: 


Copper 

0.25 

Iron: 

0.04 

Manganese: 

0.35 

Carbon: 

0.15 

Silicon: 

0.35 

Sulfur. 

0.01 

The  structure  of  nickel  is  face-centered  cubic.  Commercial  uses  for  Nickel  200  include  the 
manufacture  of  a  variety  of  processing  equipment,  particularly  to  maintain  product  purity  in 
handling  foods,  synthetic  fibers,  and  alkalies. 

The  final  samples  should  have  a  homogeneous  distribution  of  grain  sizes,  ideally  with  an 
isotropic  configuration  to  satisfy  the  arguments  used  in  the  theoretical  development  in  section  V.C. 
The  processes  of  cold  working,  recrystallization,  and  grain  growth  were  used  to  achieve  the 
desired  product. 

Cold  work  was  achieved  in  this  application  by  rolling  the  nickel.  The  processing  steps  for 
the  nickel  must  take  into  account  the  three  steps  ctf  annealing:  recovery,  recrystallization,  and  grain 
growth.  For  the  desired  metal  samples,  the  recovery  stage  should  be  complete  and  the  process  erf 
reciystaftization  should  begin,  and  the  grain  growth  must  be  limited  in  order  that  the  distribution  of 
the  grain  size  will  contain  very  small  as  well  as  large  grains.  Six  factors  influence  recrystallization: 
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amount  of  prior  deformation,  temperature,  time,  initial  grain  size,  composition,  and  amount  of 
recovery  prior  to  the  start  of  recrystallization  [Dieter  1961  p.  155];  therefore  different  annealing 
durations  will  be  tested  to  achieve  the  desired  structure. 

There  were  three  additional  steps  used  for  the  sample  preparation:  1)  the  final  samples  used 
for  ultrasonic  attenuation  measurement  must  be  less  than  .10  inch  thick,  flat  and  parallel  to  within  ± 

^•mil,  2)  photomicrographs  of  the  grain  boundaries  are  required,  and  3)  the  final  sample  must  be 

polished  on  one  side  to  a  mirror  finish. 

All  processing  steps  were  completed  in  the  Metallographic  Laboratory  of  the  Materials 
Engineering  Department  of  Drexel  University.  The  original  nickel  plate  was  approximately  6  x  5  x 

finches,  and  was  initially  annealed  and  rolled  by  Inco  Alloys.  An  analysis  was  completed  on  the 

nickel,  using  polishing  and  etching  to  bring  out  the  grain  boundaries.  The  photomicrographs 
showed  that  the  grains  were  deformed  from  Inco’s  processing  steps.  It  was  therefore  determined 
that  an  additional  rolling  and  annealing  process  was  necessary  to  allow  recrystallization  to  occur, 
resulting  in  the  grain  shape  and  size  distributions  desired.  The  nickel  was  cut  into  1-inch  wide 
strips,  and  then  the  rolling  was  done  on  a  2  roll  Stanat  mill.  Unidirectional  rolling  was  used  while 
the  nickel  strip  was  fed  in  for  each  pass,  reversing  sides  each  time  to  minimize  warpage.  Rolling 
reduction  per  pass  were  approximately  10  mils,  accumulating  in  a  final  reduction  of  50%. 

The  rolled  strip  was  cut  into  approximately  lxl  inch  tiles  on  a  slow  speed  band  saw.  The 
samples  were  then  mechanically  ground  to  within  ±1/2  mil  tolerance  for  parallel  tile  faces. 

The  tiles  were  first  used  to  establish  the  annealing  times  needed  to  cause  recrystallization  to 
occur.  One  sample  was  annealed  at  800  °C  for  1  minute;  the  steps  for  the  anneal  began  with 
opening  the  furnace  door,  and  the  samples  were  placed  on  a  steel  tray  which  was  at  temperature 
(800  °C).  Dwell  times  (in  this  case,  1  minute)  began  once  furnace  door  closed.  The  thermocouple 
reading  of  furnace  atmosphere  dropped  to  about  720  -  730  °C  while  the  furnace  door  was  open. 
The  steel  tray  was  used  to  minimize  heat  loss  upon  opening  oven  doors.  The  samples  were  air 
cooled  on  insulating,  room  temperature  ceramic  blocks. 

After  polishing,  etching  50/50  nitric/glacial  acetic  acid  mixture,  and  making  a  micrograph,  it 
was  clear  that  the  recovery  stage  of  annealing  had  not  been  completed  and  recrystallization  had  not 
begun.  The  other  two  tiles  were  annealed  for  2  and  5  minutes,  and  the  micrographs  showed  that 
the  grain  sizes  desired  is  achieved  for  these  annealing  durations.  Thus  the  three  remaining  tiles  not 
yet  processed  were  annealed  for  2, 3  and  5  minutes. 

The  next  step  was  to  obtain  high  quality  micrographs  of  the  grain  boundaries  of  the  final 
three  tiles.  To  achieve  this,  the  samples  were  ground  on  120  -  280  -  320  -  400  -  600  grit  papers  in 
order  to  remove  previous  grit  paper's  deformation  layers.  They  were  then  polished  on  a  canvas 
cloth  impregnated  with  9  micron  diamond  slurry,  followed  by  a  0.3  micron  alumina  polish  on  a 
Selvyt  cloth. 

To  bring  out  the  grain  boundaries,  the  samples  were  immersion  etched  for  40-45  seconds 
in  a  50/50  nitric/glacial  acetic  acid  mixture.  The  alumina  polish-etch  cycle  was  repeated  four  more 
times  to  remove  as  much  grinding-induced  deformation  as  possible. 

After  the  final  etch,  photographs  were  taken  on  a  Zeiss  ICM  405  Metallograph.  33 
photographs  in  total  were  taken:  1 1  at  each  annealing  time,  9  at  lOOx  for  analysis  and  one  each  at 
32x  and  200x.  The  nine  micrographs  at  lOOx  were  used  for  grain-size  distribution  measurement  at 
the  Lewis  Research  Center  of  the  National  Aeronautics  and  Space  Administration  (NASA);  the 
image  processing  technique  is  described  in  section  II.  A  final  alumina  polishing  was  done  to 
remove  the  etched  surface  and  create  a  mirror  finish.  The  final  thicknesses  were  measured  using  a 
micrometer  to  be  0.1309,  0.1299  and  0.1279  centimeters  for  the  two-,  three-  and  five-minute 
samples,  respectively. 
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After  the  ultrasonic  measurements  were  complete,  additional  photomicrographs  were  made 
on  the  side  surfaces  of  the  sample.  These  additional  micrographs  were  taken  to  establish  if  the 
grains  are  equiaxed. 

V.C.  Attenuation  Measurement 

For  metals,  the  standard  method  of  attenuation  measurement  involves  the  transmission  erf"  a 
pulse  of  ultrasonic  energy  into  the  sample,  followed  by  the  recording  of  multiple  echoes  of  the 
pulse  as  it  travels  across  the  specimen  and  is  reflected  at  the  edges  of  the  material.  Fourier 
transforms  of  these  echoes  allow  the  computation  of  the  attenuation  for  a  range  of  frequencies. 

Attenuation  measurements  reported  here  were  made  at  the  Lewis  Research  Center  of  the 
National  Aeronautics  and  Space  Administration  (NASA).  The  equipment  used  for  these 
measurements  is  very  similar  to  the  equipment  used  for  the  preliminary  ultrasonic  experiments 
using  stainless  steel  samples,  which  were  completed  in  the  Signal  Processing  Laboratory  at  Drexel 
University.  Therefore,  it  is  possible  to  make  the  attenuation  measurements  in  the  Signal 
Processing  Laboratory.  However,  since  it  was  necessary  to  make  the  grain-size  distribution 
measurements  at  the  Lewis  Research  Center,  the  decision  was  made  to  complete  the  ultrasonic 
analysis  there  as  well,  with  a  site  visit  by  D.  Nicoletti. 

Two  major  obstacles  to  accurate  ultrasonic  attenuation  measurements  are  diffraction  and 
reflection  coefficient  effects,  as  well  as  the  need  for  the  attenuation  to  be  small  enough  to  allow 
accurate  ultrasonic  recordings.  The  experimental  design  presented  below  strikes  a  balance  among 
these  considerations. 

The  samples  are  very  thin  (less  than  .  15  centimeters  thick)  when  compared  to  metal 
samples  used  in  other  attenuation  measurements.  For  comparison,  previous  researchers  have  used 
1  cm  [Smith  et.  al.  1981],  1.9  cm  [Papadakis  1963],  0.63  to  2.9  cm  [Klinman  &  Stephenson 
1981].  A  30  MHz  center  frequency  transducer  was  used.  The  attenuation  measurement  scheme 
for  pulse-echo,  contact  mode  [Generazio  1985],  is  shown  in  Figure  3. 

A  quartz  buffer  rod  and  glycerin  couplant  are  used  between  the  piezo-electric  crystal  and  the 
specimen.  Each  echo  (FSi,  FS2,  Bj,  and  B2)  was  obtained  by  performing  64  time  averages.  The 
echoes  were  first  recorded  without  the  sample:  FSi,  from  the  front  surface  of  the  buffer  rod.  Once 
the  transducer  was  in  contact  with  the  specimen,  the  echo  was  again  recorded  from  the  front 
surface  of  the  buffer  rod,  labeled  FS2.  Thus  the  frequency  dependent  buffer-couplant-sample, 
R(f),  can  be  obtained  using  the  magnitudes  of  Fourier  transforms  of  FSi  and  FS2: 

IR(OI=nilW  <38> 

FSi(f)  and  FSgfQ  indicate  the  magnitude  spectra  of  FSi  and  FS2,  calculated  using  a  Fast  Fourier 
Transform  (FFT).  With  the  transducer  still  in  contact  with  the  specimen,  the  first  and  second 
echoes  from  within  the  specimen,  Bj  and  B2,  were  recorded.  The  magnitude  of  the  Fourier 
transforms  of  these  echoes  can  be  expressed  as 

IBi  (01  =  IFSi(0l  (1  -  IR(0l2)  exp  {-2  c<0  d> 

©2  (01  =  IFSi(0l  (1  -  IR(0l2)  IR(0lexp  {-4  o(0  <0 
where  d  is  the  thickness  of  the  sample. 

The  attenuation  can  be  found  from  these  four  measurements: 


(39) 


(40) 
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where  ln()  denotes  the  natural  logarithm  function.  The  resulting  units  few  the  attenuation  is  Nepers 
per  centimeter  (Np/cm). 

Piezoelectric  transducers  are  in  general  inherently  band-limited  devices,  therefore  the  usable 
range  of  frequencies  (or  wavelengths)  for  the  measured  attenuation  must  be  chosen.  The 
frequency  range  for  the  above  measured  attenuation  is  a  result  of  the  FFT  calculations,  and  is  a  set 
of  discrete  points  from  0  Hz  to  half  the  sampling  rate  (2.56  GHz).  An  upper  limit  for  the 
attenuation  should  be  set  to  assure  that  the  measured  echoes  are  not  attenuated  to  the  extent  that 
their  measurement  is  overpowered  by  the  noise  level.  An  upper  limit  of  5  Np/cm  was  chosen.  For 
this  value  of  attenuation,  given  a  sample  thicknesses  of  0. 12  cm,  the  amplitude  of  the  second  echo 
is  e-5*3*12  or  0.016  times  the  input  amplitude.  Since  attenuation  is  inversely  related  to  wavelength 
(Equation  9),  an  upper  limit  of  5  Np/cm  gives  the  lower  limit  for  the  wavelength  range.  Other 
criteria  can  be  used  for  the  choice  of  the  usable  range  of  wavelengths,  including  the  3  dB  or  6  dB 
bandwidth  of  the  transducer,  or  consideration  of  the  experimental  uncertainty  associated  with  high 
attenuation  ultrasonic  measurement 

The  attenuation  theoretically  should  decrease  or  remain  constant  with  increasing 
wavelength.  The  upper  limit  for  the  wavelength,  X^,  will  be  chosen  to  include  the  attenuation 
curve  up  to  the  point  it  starts  to  increase  with  wavelength;  i.e.,  if  a(Xi+i)  >  a(Aj)  then  i  =  N  and 
a(Xj)  is  the  last  value  used  in  the  presented  attenuation  measurement.  The  upper  and  lower  limits 
are  expected  to  be  change  from  sample  to  sample.  However,  in  order  to  compare  the  attenuation 
curves  for  all  three  samples,  the  values  for  the  lower  and  upper  limit  for  the  wavelengths  are 
chosen  using  the  attenuation  measured  for  the  two-minute  sample. 


VI.  Experimental  Results 

A.  Grain-Size  Distribution 

1.  Experimental  Results 

The  experimental  results  of  the  grain-size  distribution  measurement  are  presented  in  this 
section.  The  isotropy  of  the  distribution  is  investigated,  and  the  mean  and  variance  of  the 
distributions  are  calculated.  The  inverse  power-law  approximation  is  displayed  and  statistically 
analyzed. 

The  log-normal  distribution  has  been  previously  used  to  describe  the  grain-size  distribution 
in  three-dimensions.  It  is  worth  repeating  that  a  log-normal  distribution  in  three-dimensions  does 
not  imply  that  the  distribution  in  two-dimensions  is  also  log-normal.  The  two-dimensional 
representation  of  the  grains  from  the  micrographs  are  used  in  this  work  for  the  size  distribution 
measurement  Therefore,  it  is  inappropriate  to  compare  the  log-normal  fit  of  the  data  to  the  inverse 
power-law  fit  since  testing  the  appropriateness  of  the  log-normal  fit  in  two-dimensions  has  no 
bearing  on  a  possible  log-normal  fit  of  the  three-dimensional  grain-size  distribution. 

The  grain-size  distribution  was  measured  at  the  NASA  Lewis  Research  Center  from  nine 
photomicrographs  of  each  sample,  using  the  tone-pulse-encoded  technique  described  in  section  II. 
Using  the  notation  from  section  IV.C,  the  histograms  are  a  measurement  of  N^(d),  defined  as 

N(d)  dd  =  number  of  grains  of  size  between  d  and  d  +  dd 

These  histograms  are  shown  for  the  three  samples  in  Figure  4a.  Each  histogram  consists 
of  256  data  points;  the  67  non-zero  points  (Ng(d)  *  0.0)  are  used  for  the  figure.  As  discussed  in 
section  II,  the  algorithm  used  for  the  grain-size  distribution  measurement  has  a  lower  limit  for  the 
resolution  of  small  grains  due  to  the  need  for  at  least  ten  [axels  to  encode  a  three-cycle  sine  wave. 
This  lower  limit  explains  the  abrupt  cut-oIT  of  the  grain-size  distributions  for  small  grain  sizes. 
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The  average  grain  size  can  be  computed  from  the  histogram  as  follows: 


average  =  E[d] 


and  are  as  follows: 


255 

l  diNaOi) 

M) _ 

“=  255 

l  NaOi) 

i=0 


Two-minute  sample: 
Three-minute  sample: 
Five-minute  sample: 


Average 
88.0  microns 
90.8  microns 
96.7  microns 


(41) 


As  expected,  the  average  grain  size  grows  with  increased  annealing  time,  also  seen  by  a 
relative  comparison  of  the  histograms  in  Figure  4a. 

The  grain-size  distributions  and  averages  are  very  close  for  the  three  samples  -  the  average 
grain  sizes  differ  by  as  little  as  3.2  per  cent  This  is  desirable  due  to  the  wavelength  considerations 
in  the  attenuation  measurement.  In  order  to  measure  the  grain  size  exponent,  the  wavelength  used 
must  be  within  the  range  of  the  grain  sizes.  Also,  in  general  transducers  inherently  have  finite 
bandwidths.  Therefore,  if  samples  were  prepared  with  very  different  grain  size  ranges,  it  would 
not  be  possible  to  use  a  single  transducer  for  the  ultrasonic  measurements.  The  alternative  is  to  use 
multiple  transducers  with  different  center  frequencies,  introducing  the  problem  of  comparing 
measurements  from  different  transducers. 

VI.A.2.  Inverse  Power-Law  Approximation 

The  results  of  the  experiment  designed  to  test  the  hypothesis  presented  in  section  IV.B  are 
presented  below,  namely  whether  a  power-law  dependence  is  a  suitable  model  for  the  grain-size 
distribution.  The  histograms  are  presented  on  a  log-log  scale  in  Figure  4b.  Except  for  the  small 
number  (five)  of  assumably  anomalous  points  at  the  beginning  and  end  of  the  histograms,  a  the 
power-law  approximation  seems  appropriate. 

There  are  two  possible  explanations  for  these  five  points  in  question.  As  discussed  in 
section  IV.C,  the  inverse  power-law  can  be  valid  over  a  finite  range  of  sizes  only.  Therefore,  one 
explanation  for  these  five  points  is  that  they  lie  outside  of  the  region  described  by  the  inverse 
power-law,  [Dp  D2],  For  this  case,  the  exponent  for  the  wavelength  dependence  of  the 
attenuation  is  determined  by  the  grain-size  distribution  exponent  as  long  as  the  inverse  power-law 
dominates  the  overall  size  distribution.  A  second  possible  explanation  comes  from  a  property  of 
the  algorithm  used  to  calculate  the  size  distribution.  As  mentioned  in  section  II,  the  grain  size 
measurement  for  the  smaller  grains  is  expected  to  have  a  high  degree  of  uncertainty  than  the 
measurement  for  the  larger  grains.  Thus,  the  difference  between  the  measurement  for  the  three 
points  at  the  beginning  of  the  curves  and  the  power-law  approximation  may  be  due  to  this 
algorithm-based  error. 

The  mean-square-error  (MSE)  and  correlation  coefficient  are  used  for  a  quantitative 
judgement  These  statistical  measurements  were  applied  to  the  entire  data  set  For  comparison, 
calculations  were  also  made  for  the  set  excluding  the  three  points  at  the  beginning  and  the  two 
points  at  the  end  of  the  curve;  this  subset  of  the  data  points  is  labeled  “partial  data  set”  below. 

The  grain-size  distributions  measurements  with  the  lines  indicating  the  linear 
approximations  based  on  the  entire  data  set  are  shown  in  Figure  5a,  6a  and  7a  for  the  two-,  three- 
and  five-minute  samples,  respectively.  Linear  regression  was  used  to  calculated  the  slope  and  y- 


19 


intercepts.  For  comparison,  the  linear  approximations  based  on  the  partial  data  set  are  presented  in 
Figure  5b,  6b  and  7b  for  the  two-,  three-  and  five-minute  samples,  respectively. 


Entire  data  set: 

Grain  size  range:  101 29  to  10244  (19.5  to  ITT)  microns 


slope 

v-intercept 

Dower-law  expression 

Two-minute  sample: 

-1.24 

4.15 

Na(d)  =  14,200  a124 

Three-minute  sample: 

-1.18 

4.07 

Na(d)  =  11,700  d’1'18 

Five-minute  sample: 

-1.04 

3.88 

Na(d)  =  7,640  a'1 04 

Partial  data  set: 

Grain  size  range:  101 49  to  102 43 

(31.2  to  269)  microns 

slope 

v-interceDt 

power-law  expression 

Two-minute  sample: 

-1.45 

4.63 

Na(d)=  42,900  a*145 

Three-minute  sample: 

-1.38 

4.54 

Na(d)  =  34,600  d1 38 

Five-minute  sample: 

-1.22 

4.30 

Na(d)  =  20,100  aL22 

To  measure  the  appropriateness  of  the  linear  relationship  for  the  log-log  representation 
the  attenuation  data,  the  sample  correlation  coefficient  r  was  computed: 

Entire  data  set: 

Grain  size  range:  101 29  to  102 44 

( 19.5  to  277)  microns 

r 

E? 

Two-minute  sample: 

-0.925 

85.6% 

Three-minute  sample: 

-0.918 

84.3% 

Five-minute  sample: 

-0.894 

80.0% 

Partial  data  set: 

Grain  size  range:  101 49  to  102  43 

(31.2  to  269)  microns 

r 

r2 

Two-minute  sample: 

-0.994 

98.7% 

Three-minute  sample: 

-0.989 

97.9% 

Five-minute  sample: 

-0.988 

97.5% 

The  negative  sign  for  the  sample  correlation  coefficient  is  due  to  the  negative  slope.  This  statistical 
measurement  shows  that  over  97%  of  the  variations  in  the  values  of  log(distribution)  are  accounted 
for  by  a  linear  dependence  on  log(size)  for  the  partial  data  set 
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A  second  method  of  measuring  the  closeness  of  the  linear  fit  is  through  the  mean  square 
error  (MSE)  and  root  mean  square  error  (RMSE).  The  results  are: 

Entire  data  set: 


Grain  size  range:  101 29  to  102  44 

( 19.5  to  277)  microns 

MSE 

RMSE 

Two-minute  sample: 

0.00157 

0.0397 

Three-minute  sample: 

0.00156 

0.0505 

Five-minute  sample: 

0.00235 

0.0485 

Partial  data  set: 

Grain  size  range:  101 49  to  102  43 

(31.2  to  269)  microns 

MSE 

RMSE 

Two-minute  sample: 

0.00170 

0.0412 

Three-minute  sample: 

0.00255 

0.0505 

Five-minute  sample: 

0.00235 

0.0485 

As  outlined  in  section  IV.A,  self-similar  distributions  are  invariant  to  magnification.  In  an 
attempt  to  demonstrate  this  for  the  grain-size  distributions,  the  photomicrographs  were  magnified 
using  a  photocopier.  The  results  were  unsatisfactory  because  the  distribution  is  approximately  1.2 
decades  wide,  a  range  not  broad  enough  to  allow  self-similarity  analysis  over  many 
magnifications.  Therefore  the  appropriateness  of  the  inverse  power-law  for  the  grain-size 
distribution  was  judged  using  the  measured  histograms  alone. 


VI.B.  Attenuation  Measurement 

1.  Experimental  Results 

The  attenuation  experiments  were  carried  out  at  the  NASA  Lewis  Research  Center. 
Measurements  were  made  from  different  transducer  locations  for  each  sample:  six  locations  for  the 
two-minute,  three  locations  for  the  three-minute,  ami  four  locations  few  the  Five-minute  samples. 
The  choice  for  the  number  of  locations  for  each  samples  was  arbitrary;  more  measurements  were 
made  for  the  two-minute  sample  since  it  was  the  first  sample  used  for  the  ultrasound  attenuation 
experiment.  The  averaged  44-point  curves  are  shown  in  Figure  8  for  the  three  samples,  for  both 
the  linear-linear  and  log-log  axes.  The  curves  are  approximately  linear  on  the  log-log  scale, 
indicating  that  a  power-law  dependence  on  wavelength  is  an  appropriate  model  for  attenuation. 

Before  the  attenuation  exponent  is  presented,  it  is  interesting  to  note  in  Figure  8b  that  there 
is  a  “knee”  in  the  data  presented  at  approximately  102  4  or  250  microns.  This  observation  is 
explained  by  comparing  the  grain  size  and  wavelength  ranges.  The  limits  for  grain  size  are  101-3 
and  102-4  (20  and  250)  microns,  while  the  wavelength  range  from  1021  to  102  5  (125  -  310) 
microns.  As  discussed  in  section  IV.C,  the  wavelength  must  be  between  the  largest  and  smallest 
grain  size  in  order  to  measure  the  size  distribution’s  exponent  Thus,  the  wavelengths  are  equal  to 
or  greater  than  the  largest  grain  size  for  the  region  \  =  lO2  4  to  102  5  (250  to  3 10)  microns,  and  the 
change  in  slope  is  believed  to  be  simply  the  result  erf'  this  fact  As  the  wavelength  becomes  much 
greater  than  the  largest  grain  size,  a  crossover  occurs  and  the  exponent  should  approach  4,  the 
theoretical  limit  for  Rayleigh  scanning. 

One  further  piece  of  information  that  can  be  extracted  from  the  attenuation  curves  is  the 
following:  for  which  sample  is  the  comer  more  pronounced?  This  is  answered  theoretically  by 
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deciding  which  sample  is  dominated  more  by  Rayleigh  scattering  for  the  upper  region  of  the  curve, 
i.e.  for  the  wavelength  range  of  102  4  to  102  5  (250  to  310)  microns.  Because  of  the  relative 
shapes  for  the  histograms  (Figure  4),  the  wavelengths  for  this  region  will  be  greater  than  a  larger 
percentage  of  the  entire  number  of  grains  for  the  two-minute  sample  than  it  is  for  the  three-  and 
five-minute  samples.  It  is,  therefore,  anticipated  that  the  comer  will  be  more  pronounced  for  the 
two-minute  sample;  this  is  also  exhibited  in  Figure  8. 

VI.B.2.  Inverse  Power-Law  Approximation 

The  exponents  were  computed  for  the  wavelength  range  up  to  102  4  or  250  microns  using 
linear  regression  to  find  the  slope.  The  attenuation  measurements  with  the  lines  indicating  the 
linear  approximations  are  shown  in  Figures  9,  10,  and  11  for  the  two-,  three-,  and  five-minute 
samples,  respectively.  A  line  with  a  slope  of  -4  is  drawn  for  the  range  102  4  to  102  5  (250  to  310) 
microns  to  show  that  the  exponent  is  approaching  the  Rayleigh  scattering  exponent  (4). 

Below  is  a  table  for  the  computed  slopes  and  y-intercepts  for  the  data  points  from  the 
wavelength  region  up  to  102  4  or  250  microns: 


slope 

v-intercept 

power-law  expression 

Two  minute  sample: 

-2.401 

5.6937 

a(X)  =  493,969  X'2  401 

Three  minute  sample: 

-2.339 

5.6253 

a(X)  =  421,988  X'2339 

Five  minute  sample: 

-2.338 

5.7078 

a(X)  =  510,270  X'2338 

These  values  for  the  exponent  are  clearly  within  the  range  previously  measured  for  Nickel  200 
samples  using  very  similar  preparation:  1.6  to  2.6  [Generazio  1988c].  The  exponents  decrease 
(slightly)  as  the  annealing  duration  increases,  also  seen  in  the  results  reported  in  Generazio  1968c. 

Error  bars  are  provided  in  Figures  9, 10,  and  1 1  for  selected  data  points.  The  experimental 

error  was  computed  ns  follows.  The  standard  deviation  (a)  was  found  for  the  linear-linear 
representations  of  the  attenuation  curves  for  each  wavelength.  Therefore  on  the  linear-linear  scale, 

the  average  attenuation  has  an  experimental  uncertainty  of  ±  o.  To  represent  the  error  on  the  log- 
log  scale,  the  error  bar  was  defined  as  follows.  For  a  given  attenuation  value  of  afX),  the  error  bar 

was  computed  as  either  |log(a(X)  +  a)  -  log  (a  w)}  or  {log(a(X)-a)-log(  a(X))j ,  whichever 
was  larger  in  magnitude. 

To  measure  the  appropriateness  of  the  linear  relationship  for  the  log-log  representation  of 
the  attenuation  data,  the  sample  correlation  coefficient  r  was  computed  using  the  wavelength  range 
up  to  102  4  or  250  microns: 

l  C2 

Two  minute  sample:  -0.9989  99.77% 

Three  minute  sample:  -0.9987  99.75% 

Five  minute  sample:  -0.9988  99.76% 

The  negative  sign  for  the  sample  correlation  coefficient  is  due  to  the  negative  slope.  This  statistical 
measurement  shows  that  over  99.7%  of  the  variations  in  the  values  of  iog(attenuation)  are 
accounted  for  by  a  linear  dependence  on  log(wavelength). 
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Similar  to  the  analysis  for  the  grain-size  distributions,  second  method  of  measuring  the 
closeness  of  the  linear  fit  of  the  log-log  attenuation  curves  is  through  the  mean  square  error  (MSE) 
and  root  mean  square  error  (RMSE);  the  results  are: 


MSE 

RMSE 

Two-minute  sample: 

0.0001021 

0.01010 

Three-minute  sample: 

0.0001053 

0.01026 

Five-minute  sample: 

0.0001024 

0.01012 

VI.C.  Comparison  of  Grain  Size  and  Attenuation  Exponents 

A  direct  superposition  of  the  log-log  representations  of  the  grain-size  distribution 
measurements  and  the  attenuation  measurements  to  test  this  hypothesis  is  not  possible  for  two 
reasons.  First,  as  mentioned  in  section  I  V.D  the  ultrasonic  attenuation  measurement  is  performed 
in  three  dimensions,  while  grain-size  distribution  is  extracted  from  a  two-dimensional 
photomicrograph.  The  exponent  for  the  attenuation  must  be  reduced  to  the  lower  dimension  by 
subtracting  one  before  comparing  it  to  the  grain-size  distribution.  This  correction  factor  could  be 
made  using  the  following  logic:  if  a  line  is  described  as 

y  =  mx  +  b  (42) 

then  a  second  line,  y',  with  a  slope  of  (m  -  1)  can  be  drawn  using 

y' =  (m-l)x  +  b  =  mx  +  b-x  =  y-x  (43) 

Thus,  the  attenuation  measurement,  a(X),  could  be  corrected  by  plotting  (log{a(A.)}  -  log}.) 
versus  logX;  the  slope  of  this  curve  is  one  less  than  the  slope  of  log{a(X.)}  versus  logX.  However, 
this  leads  to  the  second  reason  for  not  juxtaposing  the  two  measurements:  the  ranges  of  wavelength 
and  grain  size  are  not  similar  enough  to  make  the  graphs  meaningful. 

Since  the  ultrasonic  measurement  was  intended  as  a  nondestructive  method  of  measuring 
the  grain-size  distribution,  the  exponent  for  the  attenuation  was  computed  for  the  wavelength  range 
of  1021  to  102-4  (125  to  250)  microns,  reduced  by  one,  and  then  lines  with  this  slope  were 
superimposed  on  the  log-  log  representations  of  the  grain-size  distribution.  The  y-intercept  of  these 
lines  were  computed  using  the  slope  predicted  from  the  attenuation  measurement  and  the  histogram 
data  points,  using  both  the  entire  set  of  data  as  well  as  the  set  excluding  the  three  points  at  the 
beginning  and  the  two  points  at  the  end  of  the  curve.  The  mean  square  error  is  used  as  the  criterion 
forjudging  the  appropriateness  of  using  the  attenuation  exponent  as  a  measurement  of  slope  of  the 
grain-size  distribution;  if  it  is,  then  the  hypothesis  that  the  exponents  for  the  attenuation  and  grain- 
size  distribution  are  identical  is  validated. 

The  exponent  computed  from  the  attenuation  curve  for  the  two-minute  sample  is  -2.401.  A 
line  with  a  slope  of  -(2.401  -  1)  =  -1.401  is  shown  superimposed  on  the  grain-size  distribution  in 
Figure  12.  Similarly,  lines  with  slopes  of  -1.339  and  -1.338  are  drawn  over  the  grain-size 
distributions  for  the  three-  and  five-minute  samples  in  Figures  13  and  14,  respectively. 

The  closeness  of  the  linear  fit  is  based  on  the  mean  square  error  (MSE).  The  MSE  values 
were  computed  for  the  grain  size  histogram  using  the  slope  predicted  by  the  attenuation  exponent. 
These  statistical  measurements  are  applied  to  the  entire  data  set  For  comparison,  calculations  are 
also  made  for  the  set  excluding  the  three  points  at  the  beginning  and  the  two  points  at  the  end  of  the 
curve;  this  subset  of  the  data  points  is  labeled  “partial  data  set”  below. 


Entire  data  set: 


Grain  size  range:  101 29  to  102  44 

(19.5  to  277)  microns 

MSE 

RMSE 

Two-minute  sample: 

0.0256 

0.160 

Three-minute  sample: 

0.0256 

0.160 

Five-minute  sample: 

0.0325 

0.180 

Partial  data  set: 

Grain  size  range:  10149  to  102  43 

(31.2  to  269) 

microns 

MSE 

RMSE 

Two-minute  sample: 

0.00184 

0.0429 

Three-minute  sample: 

0.00268 

0.0517 

Five-minute  sample: 

0.00318 

0.0564 

Theoretically,  the  entire  size  distribution  should  be  linear  with  a  slope  equal  to  the 
attenuation  exponent  minus  one;  the  upper  portions  of  the  curve  in  Figure  12  and  14  appear  to  be 
better  fits.  To  demonstrate  this  quantitatively,  the  MSE  and  RMSE  calculations  were  made  as 

above  except  the  lower  limit  range  of  grain  sizes  used  in  the  calculations  is  changed  to  io2  006  or 
102  microns.  These  values  are: 

Grain  size  range:  102  01  to  10243  (102  to  269)  microns 


MSE 

RMSE 

Two-minute  sample: 

0.00167 

0.0409 

Three-minute  sample: 

0.00288 

0.0537 

Five-minute  sample: 

0.00169 

0.0411 

The  fact  that  the  upper  portion  of  the  curve  is  a  better  fit  (for  the  two-  and  five-minute 
samples)  is  not  a  discouraging  sign  for  two  reasons.  First  of  all,  the  upper  portion  contains  a 
larger  percent  of  the  total  data  points  because  the  grain  size  histogram  was  originally  generated  on 
the  linear  scale  and  then  represented  on  the  logarithmic  scale.  Secondly,  as  discussed  in  section  II, 
the  algorithm  used  to  measure  the  size  distribution  is  more  accurate  for  the  larger  grain  sizes 
[Generazio  1988b].  Therefore,  the  hypothesis  that  the  attenuation  exponent  is  a  valid  and 
nondestructive  measurement  of  the  exponent  for  the  grain-size  distribution  is  validated  with  this 
experimental  evidence. 

VI.D.  Scaling  Relationships  for  Annealing  Times 

This  section  examines  the  scaling  of  both  the  grain-size  distribution  and  ultrasonic 
attenuation  as  the  annealing  time  increases.  If  one  inspects  the  effect  of  the  annealing  time  on  the 
size  distribution  (Figure  4),  it  appears  that  the  curves  simply  shift  upwards  with  increased 
annealing  time,  and  the  shape  of  the  distribution  does  not  change  dramatically. 

This  observation  is  reinforced  by  a  metallurgical  process  known  as  normal  grain  growth 
[Atkinson  1988].  This  type  of  growth  is  observed  in  a  wide  variety  of  materials,  and  has  two 
properties:  1)  a  relatively  narrow  range  of  grain  sizes  and  shapes,  and  2)  a  time-invariant 
distribution  of  sizes.  In  other  words,  as  grain  growth  proceeds  with  increased  annealing  time,  the 
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form  of  the  distribution  is  preserved.  These  two  properties  are  exhibited  in  the  size  distributions 
measured  from  the  nickel  samples  used  in  this  study. 

This  section  assumes  that  the  size  distribution  simply  scales  as  the  annealing  time  grows. 
This  approach  was  also  examined  in  Generazio  1986,  using  an  average  grain  size  measurement 
instead  of  the  grain-size  distribution.  Thus  basis  curves  could  be  formed  from  the  known 
attenuation  and  grain-size  distribution  for  one  sample  (either  a  priori  or  through  a  destructive 
measurement),  then  the  size  distributions  for  other  samples  could  be  found  by  making  relative, 
scaling  measurements  of  the  attenuation.  The  attenuation  and  grain-size  distribution  measurements 
for  the  two-minute  sample  comprise  the  basis  curves,  and  the  three-  and  five-minute  curves  are 
approximated  below  by  making  the  appropriate  shifts. 

If  Nj(D)  is  the  histogram  for  the  i-th  minute  sample  (i=2, 3  or  5),  then  this  scaling  can  be 
represented  with  inverse  power-law  with  identical  exponents: 

N2(D)  =  Kj  D‘Y;  N3(D)  =  K3  D'y;  N5(D)  =  K5  D*  (44) 

The  slopes  of  the  log-log  representations  of  these  functions  will  in  theory  be  identical  for  all  three 
samples,  and  equal  to  -  y.  N2(D)  can  be  shifted  by  jlog^J-logfKj)}  to  form  N3(D). 

Similarly,  N5(D)  is  constructed  by  shifting  N2(D)  by  |log(Ks)  -  log(K2)j . 

It  has  been  proven  [Nicoletti  1991]  that  the  shifts  for  the  distribution  of  circles  on  the 
cutting  plane  are  identical  to  the  shifts  for  the  distribution  of  spheres  in  a  volume.  Na.2(d)  is 

shifted  by  log(K~)  -  log(K-)  to  obtain  N  -(d),  and  similarly  N  -(d)  is  the  result  of  shifting 
N^bylog^-logd^). 

Equation  18  gives  the  following  expression  for  the  attenuation: 

00 

a(X)  =  K  X~yJx~y  a(x)  dx  (45) 

with  a  grain-size  distribution  equal  to  KD~Y.  Using  the  expressions  for  the  grain-size  distribution 
in  Equation  44,  the  three  attenuation  functions  are 

00 

a^X)  =  K 2X_yJx_y  a(x)  dx 


00 

cx^X)  =  K 3\"yJx~y  a(x)  dx 
00 

a^X)  =  K5X"y  Jx'y  a(x)  dx  (46) 

for  the  two-,  three-  and  five-  minute  samples,  respectively.  Thus  the  shifts  for  the  log-log 
representations  of  the  attenuation  measurements  are  jlog^)-  log(K2)j  and 

|  log(Ks)  -  log(K2)} ,  the  same  shifts  used  for  the  grain-size  distributions.  Therefore,  if  the  size 

distribution  exponent  (and  therefore  the  attenuation  exponent)  remains  constant  for  different 
annealing  times,  then  the  measured  shift  in  the  log-log  attenuation  curve  corresponds  to  the  shift  or 
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scaling  of  the  grain-size  distribution  as  well.  The  attenuation  measurements  resulted  in  the 
following  parameters: 


Two  minute  sample: 
Three  minute  sample: 
Five  minute  sample: 


exponent  (slope) 
-2.401 
-2.339 
-2.338 


log(K)  (y-intercept) 
5.6937 
5.6253 
5.7078 


The  exponents  (-2.401,  -2.339  and  -2.338)  are  very  close;  the  maximum  percent  change  is  only 
2.7%.  There  is  a  problem  though  when  the  y-intercepts  are  compared  for  the  two-  and  three- 
minute  samples.  The  three-minute  curve  is  higher  than  the  two-minute  curve,  yet  the  y-intercept 
computed  is  smaller  for  the  three-minute  curve.  This  is  a  result  of  the  lines  found  from  linear 
regression  intersecting  outside  of  the  range  of  the  measurement ,  forcing  logfK^)  to  be  larger  than 

logfKg).  Since  this  does  not  represent  the  observation  that  the  measured  three-minute  curve  is  an 
upward-shifted  version  of  the  measured  two-minute  curve,  the  following  correction  was  made. 
The  attenuation  curves  presented  in  Figure  8  were  translated  to  the  left  by  2.097,  creating  the 
following  new  curves: 


«?(X)  =  aj(X  -  2.097),  i  =  2,  3,  5  (47) 

The  slopes  computed  using  linear  regression  on  a*(X)  are  identical  to  the  slopes  from  a^X);  the 
y-intercept  will  change,  as  follows: 


exponent  (slope) 

log(K*)  (y-intercept) 

log(K*)-log(K*) 

a*(X): 

-2.401 

0.65905 

.... 

a*(X): 

-2.339 

0.72032 

0.06127 

a|(X): 

-2.338 

0.80570 

0.14665 

These  y-intercept  values  correctly  reflect  that  the  curves  shift  up  as  the  annealing  duration 
increases.  The  values  given  above  for  log(K*)  -  log(K|)  replace  the  shift  values  given  previously. 

Figure  15  shows  the  results  of  shifting  a2(X)  by  0.06127  to  estimate  a3(X),  and  similarly, 

shifting  a2(D)  by  0.14665  is  drawn  as  a  representation  of  a^X)  in  Figure  16.  It  is  clear  that  the 
shapes  of  the  curves  are  close  enough  to  allow  a  simple  shift  to  obtain  one  curve  from  another. 

The  hypothesis  that  the  shifts  for  the  attenuation  measurement  curves  are  equivalent  to  the 
shift  for  the  grain-size  distribution  curves  is  tested  next  N2(D)  is  shifted  by  0.06127  and  0. 14665 
in  Figures  17  and  18  to  estimate  N3(D)  and  N^D),  respectively. 

The  agreement  between  of  the  shifted  curves  and  the  measured  histograms  has  two 
important  ramifications:  first,  it  is  further  experimental  evidence  of  the  proposed  scaling 
relationships  between  attenuation  and  grain-size  distribution.  It  also  suggests  an  important 
nondestructive  evaluation  technique,  described  as  follows.  Assuming  that  the  exponent  is  the  same 
for  a  series  of  metal  samples  and  that  the  grain-size  distribution  is  (mown  for  one  sample  (either  a 
priori  or  through  a  destructive  measurement),  then  the  size  distributions  for  all  of  the  other  samples 
could  be  found  by  measuring  the  shifts  of  the  attenuation  curves. 
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VII.  Conclusions  and  Future  Directions 


This  research  deals  with  the  application  of  nondestructive  ultrasonic  testing  to  grain-size 
distribution  measurement  in  metals.  A  novel  scaling  relationship  between  ultrasonic  attenuation 
and  grain-size  distribution  is  developed,  demonstrating  how  the  complex  structure  of  the  grain 
boundaries  can  be  analyzed  using  the  measured  ultrasonic  attenuation.  It  has  been  pointed  out  that 
researchers  have  often  made  two  assumptions:  1)  Rayleigh  scattering  dominates,  and  2)  the  grain- 
size  distribution  can  be  represented  by  a  single  value,  the  average  grain  size.  The  result  of  these 
assumptions  has  been  the  use  of  the  following  expression  for  attenuation: 


a(f)~(D)V 


(48) 


This  relationship  has  two  serious  shortcomings  for  many  experimental  results,  both  of 
which  are  consequences  of  the  grain-size  distribution.  First,  the  assumption  that  the  attenuation  is 
dominated  by  Rayleigh  scattering  requires  that  the  wavelength  is  much  greater  than  the  scatterer 
size.  Since  the  largest  grain  size  may  be  significantly  larger  than  the  average  grain  size, 
wavelengths  chosen  based  on  the  average  grain  size  may  in  fact  be  comparable  to  the  largest  grain 
size.  Secondly,  dependence  on  the  third  power  of  the  average  grain  size  as  an  accurate  grain 
volume  relationship  may  not  be  valid  if  there  is  a  distribution  of  grain  sizes.  If  a  distribution 
exists,  then  the  third  moment  must  be  used,  which  in  general  is  not  equal  to  the  cube  of  the  first 
moment  (or  average). 


The  literature  survey  in  section  III  shows  that  the  only  expression  previously  used  for  a 
grain-size  distribution  function  is  the  log-normal.  The  first  contribution  to  an  alternative  theoretical 
development  for  attenuation  is  to  use  an  inverse  power-law  for  the  grain-size  distribution. 
Justification  for  this  representation  of  the  size  distribution  can  be  found  when  scaling  theory  is 
applied  to  micrographs  of  the  grain  boundaries,  as  demonstrated  in  section  IV,  as  well  as  the 
linearity  of  the  cumulative  frequency  of  the  grain  weight  on  the  log-log  scale  [Rhines  &  Patterson 
1982,  Doherty  1984]. 


The  evolution  of  this  theoretical  work  continues  to  unfold  by  showing  that  if  the  grain-size 
distribution  is  an  inverse  power-law,  then  the  wavelength  dependence  of  attenuation  is  also  an 
inverse  power-law,  with  the  same  exponent  as  the  grain-size  distribution.  Careful  attention  is 
given  to  the  limitations  of  the  new  scaling  theory  in  terms  of  a  practical  grain-size  distribution  with 
a  finite  upper  size  limit.  Thus,  a  novel  method  is  presented  relating  the  attenuation  exponent  to  the 
grain-size  distribution. 


The  dependence  of  the  measured  attenuation  exponent  on  the  relative  size  of  the  wavelength 
to  the  finite  grain  size  range  is  explored  in  section  IV.C.  If  the  experimentally  measured  attenuation 
exponent  is  four,  then  the  wavelength  is  interpreted  to  be  much  greater  than  the  largest  grain  size, 
and  the  expression  fen*  the  attenuation  becomes 


aOO 


K(P2*~Y-  Pl4~r) 

A 


(4-y)X 


\»D2 


(49) 


Conversely  if  X  «  Dj,  then  the  measured  attenuation  exponent  is  expected  to  be  zero: 


a(X)  «  K  Cd  (DfY  -  P2*Y)  X  «  Dl  (50) 

For  exponent  values  between  zero  and  four,  the  wavelength  must  be  between  the  smallest  and 
largest  grain  size,  and  the  attenuation  exponent  can  be  directly  related  to  the  grain-size  distribution 
exponent. 

One  of  the  most  significant  contributions  of  this  work  is  for  attenuation  measured  using  this 
range  of  the  relative  values  of  wavelength  and  grain  size  (Di  «  X  «  D2).  If  Rayleigh  scattering 
dominates,  then  the  dependence  on  the  fourth  power  of  wavelength  can  be  used.  After  choosing 


an  appropriate  expression  for  the  grain-size  distribution,  the  resulting  attenuation  can  be  evaluated. 
Thus  the  inverse  power-law  expression  for  the  grain-size  distribution  does  not  improve  the 
understanding  of  the  effect  of  the  size  distribution  on  the  attenuation  for  this  type  of  scattering. 
However,  when  the  wavelength  is  between  the  largest  and  smallest  grain,  the  inverse  power-law 
assumption  for  the  grain-size  distribution  provides  a  functional,  scaling  relationship  for  the  power- 
law  dependence  of  attenuation  on  wavelength.  The  hypothesis  that  a  single  exponent  can  describe 
the  grain-size  distribution  and  the  attenuation  has  not  been  proposed  prior  to  this  study,  and  the 
experimental  evidence  contained  in  this  report  serves  to  justify  this  theory. 

The  analysis  of  the  log-log  representations  of  the  histograms  indicates  that  an  inverse 
power-law  approximation  is  appropriate  for  nickel  samples.  The  range  of  sizes  measured  (19.5  to 
276  microns)  gives  the  limits  for  the  scaling  relationship  of  the  distribution. 

The  distinction  between  a  self-similar  distribution  and  an  inverse  power-law  expression  for 
the  distribution  is  an  important  point  to  canvass.  A  power-law  can  be  used  for  a  self-similar 
distribution  with  an  infinite  range  of  sizes.  It  can  also  be  used  for  a  distribution  with  finite  limits; 
i.e.,  a  distribution  which  is  not  self-similar  over  all  scales.  The  nickel  samples  used  in  this  study, 
with  a  range  of  approximately  1.2  decades,  fall  into  the  second  category  for  inverse  power-law  size 
distributions.  This  actuality  is  illustrated  by  no  observations  of  self-similarity  property  after 
examination  of  different  magnifications  of  the  micrographs,  as  discussed  in  section  VI.  A.  Another 
recommendation  for  future  work  is  to  investigate  processing  techniques  or  metals  other  than  nickel 
that  provide  size  distributions  broad  enough  to  allow  examination  of  the  grain  structure  at  many 
magnifications,  in  order  to  reveal  the  self-similarity  property  of  the  size  distribution. 

The  above  discussion  leads  to  another  aspect  of  distributions.  An  inverse  power-law  is  one 
appropriate  function  to  use  for  distributions  determined  to  be  statistically  self-similar.  If  the  size 
distribution  exists  over  many  length  scales  but  is  not  purely  self-similar  over  all  scales  as  in  the 
case  of  multifractals  [Feder  1988  p.  66],  then  other  functions  can  be  explored,  including  the  log¬ 
normal.  Section  II  reviewed  the  mathematical  connections  between  the  log-normal  and  inverse 
power-law  functions.  These  two  functions  can  be  viewed  as  two  members  of  a  class  of  functions 
describing  distributions  which  exist  over  a  multitude  of  length  scales.  The  ramifications  of  using 
other  examples  from  this  class  of  functions  or  combinations  of  power-laws  with  different 
exponents  remains  to  be  explored. 

Rather  than  forcing  a  particular  distribution  function  to  fit  a  data  set,  a  fundamental 
approach  for  choosing  one  distribution  over  another  can  be  explored  further.  A  deeper 
understanding  of  the  underlying  growth  processes  can  be  investigated  as  a  means  of  anticipating 
the  resulting  grain-size  distribution.  Also,  if  the  basic  assumptions  for  selecting  a  particular 
distribution  are  known,  then  modifications  are  possible  for  different  materials. 

A  current  research  topic  in  the  area  of  growth  processes  is  the  study  of  cellular  patterns  and 
froth.  Grain  growth  can  be  modeled  in  terms  of  the  evolution  of  soap  bubbles:  as  time 
progresses,  some  bubbles  grow  while  others  shrink,  and  the  average  number  of  sides  and  area  per 
bubble  are  functions  of  time.  This  is  an  active  topic  of  research  [Flam  1989,  Glazier  &  Stavans 
1989,  Weaire  1989],  based  on  von  Neumann’s  law  and  other  equations  relating  pressure-driven 
diffusion  in  froths.  These  pursuits  are  directly  applicable  to  the  growth  of  grains  in  metals  [Rhines 
&  Patterson  1982,  Doherty  1984].  Most  of  the  existing  analysis  of  froth  is  in  terms  of  the  time 
dependency  of  the  number  of  sides  and  area  per  bubble,  rather  than  the  size  distribution. 
Therefore,  the  connection  between  this  type  of  endeavor  and  the  results  reported  here  are  not 
straightforward  at  this  point  in  time;  nevertheless,  it  is  a  potential  area  for  further  development 

While  the  experiments  reported  here  were  designed  to  test  the  agreement  between  the 
attenuation  and  grain-size  distribution  exponents,  it  was  unanticipated  that  the  attenuation 
exponents  for  the  three  samples  would  be  nearly  the  same.  It  has  been  shown  in  a  previous  study 
that  the  attenuation  exponent  decreases  as  the  percentage  of  recrystallization  increases  [Generazio 
1988c].  It  appears  that  fix'  the  samples  used  the  extent  of  recrystallization  is  nearly  the  same  for  the 
three  samples,  resulting  in  approximately  identical  attenuation  exponents.  This  observation  allows 
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for  examination  of  possible  scaling  relationships  among  the  three  annealing  times  for  both 
attenuation  and  grain-size  distribution,  as  shown  in  section  VI. D.  One  promising  direction  is  to 
explore  the  nondestructive  evaluation  of  grain-size  distribution  using  this  property  for  a  series  of 
metal  samples  with  the  same  exponent.  If  the  size  distribution  is  known  for  one  sample  (either  a 
priori  or  through  a  destructive  measurement),  then  the  size  distributions  for  all  of  the  other  samples 
could  be  determined  by  measuring  the  relative  shifts  of  the  attenuation. 

One  possible  justification  for  assuming  that  the  exponents  does  not  vary  for  the  three 
samples  comes  from  small  angle  scattering  analysis.  It  has  been  hypothesized  that  the  exponent 
derived  from  the  small  angle  scattering  curves  is  a  signature  of  the  growth  mechanism  for  silica 
particle  aggregation,  and  is  dependent  upon  the  growth  process  [Kjems  et.  al.  1986].  While  this 
consideration  is  beyond  the  scope  of  this  report,  a  potential  future  research  direction  is  to 
investigate  the  possibility  that  the  grain-size  distribution  exponent  is  constant  for  specific  types  of 
grain  growth. 

One  of  the  drawbacks  to  the  analysis  of  the  effect  of  the  grain-size  distribution  on 
mechanical  properties  has  been  the  measurement  of  grain  size.  As  mentioned  in  section  II,  “The 
near  absence  of  recognition  of  these  features  and  of  their  effects  on  properties  is  understandable  in 
view  of  the  great  amount  of  labor  required  in  the  experimental  measurement  of  grain  volume 
distribution...”  [Rhines  &  Patterson  1982].  It  is  therefore  suggested  that  ultrasonic  attenuation 
measurement  may  prove  to  be  an  appropriate  grain-size  distribution  measurement  for  materials,  as 
it  is  shown  to  be  for  pure  nickel,  and  can  be  used  to  gain  a  deeper  understanding  of  the  grain-size 
distribution  by  physical  metallurgists. 

The  theoretical  and  experimental  work  presented  here  are  directed  toward  nondestructive 
measurement  of  the  exponent  of  an  inverse  power-law  expression  for  the  grain-size  distribution.  It 
can  be  suggested  that  die  corners  or  crossovers  of  the  attenuation,  as  discussed  in  section  V.B,  are 
approximate  measurements  of  the  upper  and  lower  size  limits.  To  accomplish  this,  a  broad  span  of 
wavelengths  is  required,  with  a  span  much  wider  than  the  grain-size  distribution  in  carder  to 
determine  the  comers.  While  this  is  experimentally  feasible,  the  problem  of  finding  the  exact 
values  for  the  limits  remains.  In  any  event,  if  the  limits  are  known,  then  any  statistical  moment  can 
be  calculated,  including  the  average  grain  size: 


D 


(53) 


Thus  the  nondestructive  measurement  of  the  largest  and  smallest  grain  size  using  very  broadband 
attenuation  measurements  is  a  future  application  of  this  theory. 

A  major  thrust  of  nondestructive  testing  has  been  flaw  detection  in  metals  using  ultrasonic 
analysis.  The  ongoing  research  implicitiy  requires  knowledge  of  grain  boundary  scattering  in 
order  to  detect  and  accurately  locate  a  flaw  site.  This  application  of  ultrasonic  analysis  treats  grain 
boundary  scattering  as  the  noise  source,  and  models  for  the  scattering  are  used  to  improve 
detection  schemes.  For  example,  the  magnitude  spectrum  for  grain  boundary  echoes  has  been 
expressed  as  [Shankar  et  al.  1986] 

Magnitude  spectrum  -  f2  e'1*  (54) 

where  b  is  a  constant.  The  derivation  of  this  expression  assumes  that  the  attenuation  due  to 
scattering  is  dependent  upon  the  fourth  power  of  frequency  and  that  the  absorption  effects  are 
significant  The  theoretical  work  continues  by  finding  an  expression  for  the  magnitude  spectrum 
from  the  flaw  signal,  and  then  band-pass  filters  and  accompanying  algorithms  are  used  to  suppress 
the  grain  noise.  More  appropriate  expressions  for  the  grain  boundary  scattering  may  improve  the 
development  and  analysis  of  the  signal  processing  techniques  used  for  flaw  detection. 

In  conclusion,  the  progress  outlined  in  this  report  makes  several  contributions  to  the 
understanding  of  attenuation  in  metals  from  a  number  of  perspectives.  The  work  presented  here 
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establishes  that  the  scaling  properties  of  grain  size  and  ultrasonic  attenuation  are  related,  and  the 
attenuation  exponent  is  shown  to  be  an  appropriate  nondestructive  measurement  of  the  grain-size 
distribution  exponent,  as  validated  by  the  experimental  evidence  using  pure  nickel.  A  further 
scaling  relationship  for  the  annealing  duration  is  shown  to  exist,  implying  the  possibility  of  more 
extensive  nondestructive  evaluation  for  samples  with  identical  exponents.  These  analysis 
techniques  have  implications  for  both  grain  size  determination  and  flaw  detection  in  metals.  Future 
endeavors  may  include  examination  of  other  candidate  materials  for  the  scaling  relationships 
offered  here;  such  an  undertaking  is  indicated  and  justified  by  the  positive  results  obtained  in  this 
study. 
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Figures 


Figure  1:  Tone-pulse  encoding  for  three  grains 


Figure  2:  Circles  of  various  sizes,  shown  at  two  magnifications 
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Figure  3:  Experimental  Set-up  for  Attenuation  Measurement 
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Number  of  grams,  linear  scale 


Figure  4a  Grain  size  histograms  for  the  two-,  three-  and 
five-minute  samples 
Linear-linear  scale 


Figure  4b  Grain  size  histograms  for  the  two-,  three-  and 
five-minute  samples 
Log-log  scale 
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Figure  5a 


Grain  size,  microns,  logarithmic  scale 

Inverse  power-law  approximation  based  on  entire 
data  set  and  measured  values  for  the  grain  size 
histogram  for  the  two-minute  sample 
Log-log  scale 


Figure  5b  Inverse  power-law  approximation  based  on 

partial  data  set  and  measured  values  for  the  grain 
size  histogram  for  the  two-minute  sample 
Log-log  scale 
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Grain  size,  microns,  logarithmic  scale 

Inverse  power-law  approximation  based  on  entire 
data  set  and  measured  values  for  the  grain  size 
histogram  for  the  three-minute  sample 
Log-log  scale 


Inverse  power-law  approximation  based  on 
partial  data  set  and  measured  values  for  the  grain 
size  histogram  for  the  three-minute  sample 
Log-log  scale 
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Figure  8a  Averaged  attenuation  measurement  for  the  two-, 
three-  and  five-minute  samples 
Linear-linear  scale 
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Figure  8b  Averaged  attenuation  measurement  for  the  two-, 
three-  and  five-minute  samples 
Log-log  scale 
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Data  from  two-minute  sample 


Figure  9  Linear  approximation  and  measured  values  for 

the  attenuation  with  emir  bars  for  the  two-minute 
sample 
Log-log  scale 


Data  from  three-minute  sample 


Wavelength,  microns,  logarithmic  scale 


Figure  10  Linear  approximation  and  measured  values  for 
the  attenuation  with  error  bars  for  die  three- 
minute  sample 
Log-log  scale 
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Attenuation,  Np/cm,  logarithmic  scale 


Data  from  five-minute  sample 


Figure  1 1  Linear  approximation  and  measured  values  for 

the  attenuation  with  error  bars  for  the  five-minute 
sample 
Log-log  scale 


Figure  12  Slope  predicted  by  attenuation  results  and 

measured  values  for  the  grain  size  histogram  for 
the  two-minute  sample 
Log-log  scale 
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Number  of  grains,  logarithmic  scale  Number  of  grains,  log.ri.hmic  scale 


Figure  13  Slope  predicted  by  attenuation  results  and 

measured  values  for  the  grain  size  histogram  for 
the  three-minute  sample 
Log-log  scale 


Figure  14  Slope  predicted  by  attenuation  results  and 

measured  values  for  the  grain  size  histogram  for 
the  live-minute  sample 
Log-log  scale 
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Number  of  grains,  logarithmic  scale  Attenuation,  Np/cm,  logarithmic  scale 


Figure  15  Measured  attenuation  for  two-minute  sample 

shifted  by  0.061 27  and  measured  values  for  the 
attenuation  for  the  three-minute  sample 
Log-log  scale 


Figure  16  Measured  attenuation  for  two-minute  sample 

shifted  by  0.14665  and  measured  values  for  the 
attenuation  for  the  five-minute  sample 
Log-log  scale 
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Number  of  grains. 


Grain  size,  microns,  logarithmic  scale 

Figure  18  Measured  grain  size  histogram  for  two-minute 

sample  shifted  by  0.14665  and  measured  values  for 
the  grain  size  histogram  for  five-minute  sample 
Log-log  scale 
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